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FOREWORD 


This  technical  note  describes  the  operational  Spectral 
Ocean  Wave  Model  (SOWM)  and  provides  results  of  several  eval¬ 
uation  studies.  The  SOWM,  operational  at  Fleet  Numerical 
Weather  Central  since  December  1974,  represents  state  of  the 
science  ocean  wave  analysis  and  forecasting  techniques.  It 
provides  an  improved  specification  of  the  sea  surface  wave 
energy  and  offers  a  challenge  to  both  the  environmentalist 
and  the  mariner  to  intelligently  use  the  increased  information 
available. 

The  SOWM  is  a  result  of  joint  efforts  of  the  Naval 
Oceanographic  Office  ana  Fleet  Numerical  Weather  Central. 

Reviewed  and  approved  15  July  1975 


Captain,  U.S.  Navy 
Commanding  Officer 
Fleet  Numerical  Weather  Central 
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ABSTRACT 


Through  the  cooperative  efforts  of  the  Naval  Oceanographic 
Office  (NAVOCEANO)  and  Fleet  Numerical  Weather  Central  (FNWC) 
personnel,  the  open  ocean  Icosahedral-Gnomonic  wave  spectral 
computer  model  developed  by  Professor  Willard  Pierson  has  been 
placed  in  real  time  operational  use  at  FNWC  under  the  title 
Spectral  Ocean  Wave  Model  (SOWM) .  The  model  computes  semi- 
daily  on-time  analyses  and  forecasts  to  seventy-two  hours. 

A  gnomonic  projection  allows  great  circles  to  be  represented  by 
straight  lines  within  the  icosahedral  triangles,  thus  minimizing 
propagation  errors. 

The  growth  of  wave  energy  for  a  given  wind  velocity  at  grid 
points  is  accomplished  by  a  modified  Miles-Phillips  techinque. 
Growth  is  limited  by  the  Pierson-Moskowitz  fully  developed  spec¬ 
trum  for  the  given  wind  speed.  The  wave  energy  direction  is 
obtained  through  an  equation  derived  by  the  Stereo  Wave  Observa¬ 
tion  Project  (SWOP) .  The  wave  energy  spectrum  at  each  grid  point 
is  represented  by  a  15  frequency  by  12  direction  matrix. 

High  wave  energy  case  studies  were  made  for  time  periods 
during  November-December  1969  and  October  1973  for  verification 
of  the  computer  model.  Measured  wave  spectra  were  obtained 
from  Ocean  Station  PAPA  (1973)  and  FLIP  (1969)  for  comparison. 

The  analyzed  wave  spectra  are  also  compared  to  measured  wave 
spectra  from  the  National  Data  Buoy  Office  (NDBO)  EB03  buoy  for 
December  1974  but  a  detailed  analysis  was  not  made  for  this 
time  period. 
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Comparison  between  the  computed  and  measured  wave  spectra 
and  analysis  of  FNWC  computed  wind  fields  indicate  that  the  SOWM 
produces  far  superior  results  than  the  previous  FNWC  operational 
"singular"  wave  model. 


AN  EVALUATION  OF  A  HEMISPHERIC  OPERATIONAL 
WAVE  SPECTRAL  MODEL 

I.  INTRODUCTION 

The  capability  of  the  U.  S.  Navy  to  analyze  and  forecast 
deep  water  wave  conditions  on  a  global  scale  has  improved 
significantly  during  the  past  year.  In  the  past,  Fleet 
Numerical  Weather  Central  (FNWC) ,  the  global  forecasting  arm 
of  the  U.  S.  Navy  Weather  Service,  has  relied  on  "singular" 
wave  models  (Hubert  and  Mendenhall,  1970;  Schwartz  and 
Hubert,  1973)  to  predict  wind  wave  and  swell  heights  and 
their  corresponding  directions  and  periods  at  points  on 
Northern  Hemispheric  grid  systems.  The  basic  weakness  of 
the  "singular"  models  is  that  they  do  not  accurately  depict 
the  complex  wave  propagation  in  the  larger  oceans  such  as  the 
Pacific  and  Atlantic  Oceans  where  several  wave  trains  can 
coexist  in  one  area  at  any  given  time. 

In  April  1972,  through  the  cooperative  efforts  of  the 
U.  S.  Naval  Oceanographic  Office  (NAVOCEANO) ,  FNWC  and  New 
York  University  (NYU  cognizant  personnel  are  now  associated 
with  the  University  Institute  of  Oceanography,  City  University 
of  New  York) ,  an  operational  Mediterranean  wave  spectral 
model  was  placed  into  use  at  FNWC  (Lazanoff,  Stevenson  and 
Cardone,  3.972)  .  The  model  computes  twice  daily  on-time 
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analyses  and  forecasts  out  to  thirty-six  hours.  Besides 
producing  a  far  superior  and  more  sophisticated  product  than 
the  "singular"  model,  the  performance  of  the  Mediterranean 
wave  spectral  computer  model  demonstrated  the  practicability 
of  running  spectral  models  on  an  operational  basis. 

It  was  then  decided  to  evaluate  the  NYU  Icosahedral- 
Gnomonic  (IG)  wave  spectral  model  as  a  real-time  operational 
hemispheric  wave  spectral  model.  The  NYU  IG  wave  spectral 
model,  modified  for  operational  use,  commenced  operational 
evaluation  in  July  1973  and  in  December  1974  was  designated 
the  operational  FNWC  hemispheric  open  ocean  wave  model  under 
the  title  Spectral  Ocean  Wave  Model  (SOWM) . 


II.  THE  ICOSAHEDRAL-GNOMONIC  WAVE  SPECTRAL  MODEL 


The  IG  wave  model  was  designed  by  Professor  Willard  J. 
Pierson  for  global  use;  however,  since  the  global  model  program 
requires  a  large  amount  of  computer  storage  and  computation 
time,  it  was  decided  that  it  would  be  impractical  to  use  the 
full  global  model  on  the  present  FNWC  computers.  The  SOWM 
is  limited  to  the  Northern  Hemisphere. 

The  origins  of  this  particular  wave  spectral  model  can  be 
traced  back  to  the  original  NYU  Atlantic  wave  model  (Baer, 

1962),  and,  as  with  the  previous  NYU  wave  models,  can  be 
divided  into  two  independent  parts:  (a)  growth  at  each  grid 
point  and  (b)  propagation  from  grid  point  to  grid  point. 

A.  Wave  Energy  Growth 

The  growth  technique  is  essentially  the  same  as  the 
one  used  in  the  Mediterranean  wave  model  and  was  originally 
developed  by  Inoue  (1967) .  Inoue  combined  the  Miles  instability 
theory  with  the  Phillips  resonance  theory  such  that  when  the 
sea  begins  to  grow  from  calm  conditions,  the  resonance  mechanism 
predominates  and  later,  as  the  wind  velocity  begins  to  increase, 
the  instability  mechanism  becomes  more  dominant.  The  Phillips 
theory  essentially  states  that  a  resonance  between  the  air-sea 
system  could  occur  when  a  component  of  the  surface  pressure 
distribution  moved  at  the  same  speed  as  a  free  surface  wave 
of  the  same  wave  number  (where  the  wave  number,  K,  is  equal 
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to  2ir/L,  and  L  is  the  wave  length) .  The  Niles  instability 
theory  states  that  the  mean  rate  of  energy  transferred  from 
the  parallel  sheer  flow  to  the  surface  wave  is  proportional 
to  the  curvature  of  the  wind  profile  at  the  height  where  the 
mean  wind  velocity  is  the  same  as  the  phase  speed  of  the  wave 
component . 

The  wave  spectra  which  are  computed  for  the  given  wind 
velocity  at  each  time  step  are  not  allowed  to  exceed  the 
Pierson-Moskowitz  fully  developed  spectrum  for  the  same  wind 
velocity.  The  directional  spectra  are  computed  by  an  equation 
developed  from  Stereo  Wave  Observation  Project  (SWOP)  data. 

The  energy  spectrum  at  each  grid  point  is  divided  into  a  15 
frequency  by  12  direction  (30°  increments)  matrix-  The 
delineation  of  the  frequency  bands  is  shown  in  Table  1.  The 
highest  frequency,  in  practice#  is  limited  from  0.164  to  0.40 
HZ  and  the  wave  energy  within  this  frequency  range  is  not 
propagated. 

One  significant  change  has  recently  been  made  to  the 
growth  portion  of  the  wave  spectral  program.  Professor  Vincent 
Cardone,  CUNY,  has  modified  the  Phillips'  resonance  mechanism 
or  initial  growth  portion  of  the  model  (personal  communication, 
1974) .  For  lower  wind  speeds  (less  than  or  equal  to  30  knots) 
the  wave  energy  will  grow  at  a  faster  rate  during  the  initial 
six  hours  using  the  Cardone  modification  than  it  did  un  \e 

old  scheme.  The  reverse  is  true  for  wind  speeds  greater  - n 

30  knots.  After  six  hours  the  new  growth  rate  is  always  slower 
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than  the  old  growth  rate.  Comparisons  between  the  two  growth 
schemes  are  shown  in  Figures  1-7  for  the  35  knot  case. 

Wave  energy  dissipation  in  this  wave  spectral  model  is  a 
function  of  the  wind -wave  interaction  (Inoue,  1967)  and  the 
land  sea  boundary  delineation.  Energy  dissipation  as  a  function 
of  wave-wave  interactions  and  breaking  waves  is  not  included; 
however,  there  are  plans  to  consider  these  parameters  at  a 
later  date.  Although  not  used  within  the  model  as  an  energy 
dissipation  mechanism,  a  whitecap  percentage  for  each  grid 
point  is  calculated.  Whitecaps  are  a  manifestation  of  wave 
breaking  and  are  considered  to  be  closely  related  to  the 
energy  transferred  to  the  fully  developed  portion  of  the  spectrum 
(Cardone,  1969) .  The  calculations  are  made  on  the  assumption 
that  the  fetch  length  is  uniform.  It  now  appears  that  the 
computation  of  whitecap  percentages  should  also  be  based  on  the 
interaction  between  wind  waves  and  swell  and  the  variability 
of  fetch  length  (Shemdin,  1973;  Ross  and  Cardone,  1974). 

Phenomena  such  as  foam  and  streaks  are  not  considered  in  the 
calculation  of  whitecap  percentages  but  may  also  prove  to  be 
useful  tools  in  accounting  for  energy  dissipation.  Knowledge 
of  the  percentage  of  whitecaps,  foam  and  streaks  in  a  geo¬ 
graphical  area  has  significance  in  other  meteorological  and 
oceanographic  studies  and  will  be  discussed  in  more  detail  in  a 
later  paper. 

Appendix  A  contains  a  list  of  the  most  significant 
equations  used  in  the  wave  spectral  model. 
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B.  Wave  Energy  Propagation 

The  most  significant  difference  between  the  previous 
wave  spectral  models  developed  by  the  Pierson  oceanography 
group  and  the  SOWM  wave  spectral  model  is  the  method  used  to 
propagate  wave  energy.  A  number  of  papers  (Lockheed,  1965, 
1966,  1967;  Pierson,  Tick  and  Baer,  1966;  Adamo,  Baer  and 
Hosmer,  1968;  Salfi,  1974,  among  others)  have  been  written, 
thoroughly  describing  the  reasons  for  selecting  this  par¬ 
ticular  grid  system  and  propagation  scheme.  Only “a  brief 
discourse  on  the  subject  will  be  presented  in  this  paper. 

Since  deep  water  wave  energy  propagates  along  great 
circles  at  group  velocity, -it  is ;highly  desirable  for  modeling, 
work  to  represent  great  circles  as  straight  lines.  The 
gnomonic  projection  is  the  only  plane  projection  which  has 
this  unique  feature.  Unfortunately,  if  this  projection  was 
to  be  used  over  large  areas  such  as  the  North  Pacific  Ocean, 
there  would  be  great  areal  distortion.  An  icosahedron  is 
used  to  depict  the  shape  of  the  globe  to  minimize  the  areal 
distortion.  An  icosahedron  is  a  twenty  sided  polygon  with 
equilateral  triangles  for  it's  faces.  Each  face  of  the  icosa¬ 
hedron  is  a  separate  gnomonic  projection.  The  icosahedral- 
gnomonic  projection  is  shown  in  Figure  8.  The  mathematical 
properties  of  the  gnomonic  projection  on  an  icosahedral  tri¬ 
angle  are  shown  in  Table  2.  Note  that  areal  distortion  is 
within  reasonable  limits  throughout  most  of  the  triangle. 

The  distortion  becomes  a  maximum  near  the  vertices.  This  is 
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demonstrated  in  Figure  9.  The  triangles  are  oriented  in  a 
manner  such  that  as  many  vertices  as  possible  are  located  on 
land;  as  many  edges  as  possible  are  along  shore  lines;  and 
as  few  faces  as  possible  are  used  for  ocean  areas. 

The  SOWM  uses  seven  triangles  for  the  North  Pacific  Ocean, 
six  triangles  for  the  North  Atlantic  Ocean,  and  one  triangle 
for  the  Indian  Ocean.  Each  triangle  has  325  grid  points 
with  a  spacing  of  approximately  350  km.  at  the  point  of  tan- 
gency  and  194  Jen  at  the  vertices.  All  points  south  of  the 
equator  are  treated  as  land  points  at  this  time.  A  time  step 
of  three  hours  is  used  in  this  model. 

The  propagation  of  the  wave  energy  is  done  by  a  modified 
velocity  gradient  technique.  For  a  given  wave  direction  along 
a  great  circle,  the  energy  as  a  function  of  frequency  is 
assumed  to  be  propagated  at  group  velocity  in  accordance  with 
linear  theory — i.e.  the  multiplication  of  the  linear  product 
of  energy  computed  from  adjacent  grid  points  by  the  group 
velocity.  It  has  been  shown  (Lockheed,  1965)  that  propagation 
by  this  linear  velocity  gradient  technique  is  inadequate  if  the 
wave  component  fields  are  not  continuous  to  the  first  order. 
Since  obvious  discontinuities  exist,  an  additional  field  is 
computed  and  updated  at  each  time  step  to  account  for  this 
problem.  Salfi  (1971)  indicates  several  criteria  under  which 
discontinuities  can  exist  in  the  energy  fields:  (a)  any  local 
maximum  or  minimum,  (b)  any  land  point  adjacent  to  a  sea  point, 


and  (c)  any  point  where  the  change  in  value  between  that  point 
and  the  down-stream  grid  point  exceeds  20%  (an  arbitrarily  set 
value)  of  the  arithmetic  average  of  the  two  points.  This 
discontinuity  field  ranges  in  value  from  zero  to  one  where 
"one"  indicates  that  no  discontinuity  exists.  Values  less  than 
one  indicate  the  fractional  distance  that  a  specified 
directional-frequency  energy  component  has  moved  towards  the 
adjacent  downstream  grid  point.  Wave  energy  in  the  highest 
frequency  band  is  considered  fully  arisen  and  is  not  propagated 
Wave  energy  in  this  frequency  band  takes  between  twelve  to 
twenty  hours  to  reach  the  adjacent  grid  point  and  thus,  it  is 
not  practical  to  propagate  this  energy.  , 

Wave  energy  as  a  function  of  direction  cannot  always  be 
propagated  directly  from  grid  point  to  grid  point  within  the 
triangles.  In  order  to  make  full  use  of  the  icosahedral- 
gnomonic  projection,  as  described  in  the  second  paragraph  of 
this  section,  the  selected  grid  system  is  based  on  a  60°  co¬ 
ordinate  system.  Remembering  that  the  wave  energy  spectra 
are  divided  into  twelve  direction  bands  and  great  circles  are 
represented  by  straight  lines  within  the  icosahedral  triangles, 
it  can  be  seen  that  two  different  propagation  situations  exist 
when  wave  energy  is  propagated  as  a  function  of  direction. 

There  are  six  primary  geometrical  directions  (o° ,  60°,  120°, 
180°,  240°,  300°)  and  six  secondary  directions  (30°,  90°,  150°, 
210°,  270°,  330°).  Wave  energy  is  propagated  directly  from 
grid  point  to  grid  point  along  the  six  primary  directions  and 
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by  zig-zag  method  along  the  six  secondary  directions.  ,  The 
zig-zag  method  is  illustrated  in  Figure  10.  Wave  energy  from 
location  3  is  propagated  on  alternate  time  steps  via  locations 
one  and  two  in  the  example.  The  same  method  is  used  for  all 
secondary  directions.  The  zig-zag  method  is  modified  near 
land  points  so  that  wave  energy  propagation  near  coast  lines 
can  be  accurately  depicted.  Near  the  coast  both  upstream 
grid  points  and  discontinuity  values  must  be  checked  at  all 
time  steps. 

A  significant  problem  arises  when  wave  energy* is  propa¬ 
gated  across  triangle  interfaces.  While  the  great  circles  are 
depicted  by  straight  lines  within  triangles,  they  refract,  as 
can  be  seen  in'  Figured,  when  crossing  from  one  triangle  to 
another  ie.,  there  is  a  discontinuity  between  triangles.  Prop¬ 
agation  across  these  interfaces  is  accomplished  by  designing 
triangles  so  that  adjacent  triangles  have  a  common  row  of  grid 
points.  When  computing  the  wave  energy  propagation  along 
this  common  row  of  points,  the  angle  of  refraction  and  the 
conversion  from  one  set  of  geometrical  coordinates  to  another 
must  be  accounted  for  because  of  the  difference  in  orientation 
between  the  triangles.  An  example  of  this  type  of  propagation 
is  shown  in  Figure  11. 

When  the  geometrical  directions  are  converted  to  geo¬ 
graphical  directions,  the  display  of  the  geographical  directions 
will  differ  from  triangle  to  triangle  for  points  located  at  the 
same  position  within  the  triangles  since  the  triangles  are  not 
all  oriented  in  the  same  direction. 
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III.  WIND  VELOCITY  INPUT 

%  — — .  ■  ■  ii  i  -i»  i 

It  cannot  be  stressed  too  strongly  that  no  matter  how 
sophisticated  and  scientifically  correct  a  wave  prediction 
model  may  be,  the  accuracy  of  the  wave  model  is  still  limited 
by  the  accuracy  of  the  wind  velocity  input  to  the  model.  The 
methods  of  calculating  wind  velocities  at  FNWC  vary  depending 
on  whether  the  program  is  in  an  analysis  or  prognosis  mode. 

At  analysis  time  the  previous  (six  hour  old)  sea  level 
pressure  (SLP)  field  is  kinematically  extracted  as  a  first 
guess.  The  SLP  fields  are  retained  on  a  hemispheric  125*125 
grid  system.  Observed  data  are  added  to  the  SLP,  by  the 
Fields  by  Information  Blending  (FIB)  technique  (Holl  and 
Mendenhall,  1971) .  The  FIB  technique  is  a  method  which  assigns 
weights  to  observed  data  depending  on  the  age  of  the  data,  the 
accuracy  of  the  observations  as  compared  to  adjacent  obser¬ 
vations  or,  in  data  spare  areas,  climatology,  and  the  source 
of  observations.  Above  20°N  and  below  20°S  fourth  order  geo- 
strophic  winds  in  u  and  v  components  are  computed  from  the 
analyzed  SLP  (Hubert  and  Mendenhall,  1970).  The  geostrophic 
wind  velocities  are  corrected  for  cuvature  and  stability  as 
a  function  of  sea  surface  temperature  advection.  A  correction 
in  proportion  to  the  sea  surface  temperature  advection  serves 
as  an  empirical  correction  for  latitude.  Persistence  (six 
hour  old  wind  velocity  analysis)  is  used  between  20°S  and 
20°N.  The  wind  analyses  are  then  modified  by  observed  wind 
data.  The  observations  must  pass  gross  error  and  lateral 
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reject  checks.  A  numerical  variation  analyses  (NVA)  scheme 
uses  the  equations  of  motion  to  couple  the  wind  and  pressure 
analyses  within  reasonable  dynamical  constraints  (lewis  and 
Grayson,  1972).  The  NVA  provides  a  method  of  extrapolating 
observations  from  data-rich  areas  to  data-def icient  areas. 

LT  Roger  Langland,  FNWC  meteorologist,  has  determined  that  the 
overall  root  mean  square  error  (RMSE)  between  the  observed 
data  and  the  final  analyses  is  on  the  order  of  2.75  meters/ 
seconds.  The  RMSE  does  appear  to  vary  depending  on  the 
amount  and  quality  of  observed  data  in  the  area  and  the  location 
of  the  area;  however,  FNWC  has  not  broken  the  errors  down 
into  finer  detail  at  the  present  time.  SLP  and  wind  field 
analyses  are- calculated  every  six  hours  on  an  operational  basis 
at  FNWC. 

The  SLP  prognoses  are  computed  by  the  FNWC  primitive 
equation  (PE)  model.  The  operational  PE  model  computes  SLP 
prognoses  every  six  hours  out  to  72  hours  on  a  twice-daily  basis. 
Modified  geostrcohic  wind  velocity  prognoses  are  then  computed 
by  the  method  described  above.  There  are  no  statistics  avail¬ 
able  for  determining  the  accuracy  of  the  wind  velocity  prognoses. 
It  is  known  that  the  PE  model  moves  typical  weather  systems  at 
approximately  82%  of  the  actual  speed  of  the  systems  because  of 
truncation  errors  in  the  second  order  differencing  scheme 
used  in  the  PE  model.  In  the  near  future  a  fourth  order 
differencing  scheme  used  in  the  PE  model  and  it  is  expected 
that  the  movement  of  the  prediction  weather  systems  will 


11 


closer  approximate  the  actual  movement  of  the  weather  systems. 

Frictional  wind  velocities  (U*)  are  computed  from  the 
analyzed  or  forecasted  wind  speeds  and  are  the  actual  input 
to  the  wave  spectral  model  growth  equations.  For  a  given 
wind  speed,  U*  can  vary  considerable  depending  on  whether 
stable,  neutral  or  unstable  conditions  exist.  An  unstable 
condition  occurs  when  the  underlying  sea  surface  temperatures 
is  warmer  than  the  above  air  temperature.  The  wave  energy 
growth  rate  is  the  highest  Under  unstable  conditions  (Cardone, 
1969) .  The  importance  of  determining  which  condition  exists 
is  shown  in  Figure  12.  As  an  example,  if  winds  were  measured 
at  20  meters,  a  wind  speed  of  about  24  knots  would  be  needed 
under  stable  conditions  to  produce  the  effective  wave  generating 
ability  of  a  wind  speed  of  only  about  19  knots  during  neutral 
conditions  and  17  knots  under  unstable  conditions  (Cardone, 

1969)  .  The  present  operational  wave  spectral  model  only 
calculates  U*  for  neutral  conditions  (Salfi,  1974) ;  thus, 
growth  rate  errors  can  exist  if  stable  or  unstable  conditions 
actually  exist.  The  greatest  errors  would  occur  under  stable 
conditions.  Note  that  the  above  example  was  based  on  wind 
measured  at  20  meters  above  the  sea  surface.  The  vertical 
change  in  wind  speed  in  the  marine  boundary  is  considered  to 
have  a  logarithmic  profile  from  the  sea  surface  to  the  top  of 
the  layer.  Stable  conditions  produce  the  more  significant 
deviations  from  the  purely  logarithmic  profile  than  equivalent 
unstable  conditions  do. 
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The  initial  growth  equation,  Phillips  Resonance  Theory, 
(Appendix  A,  Equations  3  and  4)  was  derived  from  winds  measured 
at  6.1  meters  while  the  remander  of  the  growth  equation  and 
the  Pierson-Moskowitz  fully  developed  spectra  were  derived 
from  winds  measured  at  19.5  meters.  When  Pierson  and  his  group 
use  the  wave  spectral  model,  they  put  a  great  deal  of  effort 
into  converting  the  wind  input  from  various  sources  to  the  proper 
height  for  use  in  the  growth  portion  of  the  model.  This  is  not 
presently  being  done  in  the  operational  version  of  the  IG  model. 
The  FNWC  surface  wind  fields  are  not  based  on  any  specific  height 
nor  is  any  attempt  made  to  correct  for  any  height  deviation  from 
19.5  meters  or  6.1  meters.  This  omission  may  cause  significant 
errors  under  -certain  meteorological .conditions .  The  calculation 
of  U*  based  on  stability  considerations  and  the  adjustment  of 
wind  speed  for  height  are  among  the  highest  priority  areas  of 
work  that  will  be  done  in  the  near  future  and  should  make  sig¬ 
nificant  contributions  in  improving  the  accuracy  of  the  model. 


IV.  COMPUTER  RESOURCES  REQUIRED  FOR  PRESENT  OPERATIONAL 
HEMISPHERIC  WAVE  MODEL 


The  following  CDC  6500  computer  resources  are  required  to 
run  the  SOWM  on  an  operational  basis  to  produce  an  analysis 
and  forecast  to  72  hours.  Each  time  step  requires  20  seconds 
of  Central  Processing  Unit  time  per  subprojection. 

A.  North  Pacific  Ocean  (7  Subprojections) 

1.  Computer  Time  56  minutes 

2.  Central  Memory  145000g 

3.  Extended  Core  Storage  740oK 

O 

B.  North  Atlantic  Ocean  (6  subprojections) 

1.  Computer  Time  48  minutes 

2.  Central  Memory  142000g 

3.  Extended  Core  Storage  640gK 

C.  Indian  Ocean  (1  subprojection-North  of  Equator) 

1.  Computer  Time  8  minutes 

2.  Centra'  Memory  110000g 

3.  Extended  Core  Storage  212gK 
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V.  VERIFICATION  OF  THE  WAVE  SPECTRAL  MODEL 


In  order  to  have  some  degree  of  confidence  in  any  operational 
computer  model,  the  results  of  the  model  need  to  be  verified 
with  measured  data.  Undar  normal  circumstances,  routine 
shipboard  visual  observations  are  not  adequate  for  verification 
of  a  wave  spectral  model  because  of  the  complexity  of  the  model 
output,  ship  observers’  biases,  and  the  peculiarities  of  the 
World  Meteorological  Organization  coding  system  used  to  routinely 
report,  wave  parameters.  Data  from  calibrated  wave  measuring 
instruments  such  as  Tucker  Shipboard  Wave  Recorders  and  wave 
buoys  are  considered  more  desirable.  Since  the  Icosahedral- 
Gnomonic  wave  spectral  model  was  first  used  in  an  operational- 
evaluation  mode  in  the  North  Pacific  Ocean  at  FNWC,  the  initial 
verification  studies  have  been  concentrated  in  this  area. 

Four  case  studies  are  presented.  The  first  case,  30  November  - 
2  December  1969,  demonstrates  the  ability  of  the  SOWM  wave 
spectral  model  to  propagate  low  frequency  wave  energy.  The 
second  case,  26  October  -  28  October  1973,  compares  results  of 
the  wave  spectral  model  to  measured  data  from  a  Tucker  Shipboard 
Wave  Recorder.  The  third  case,  December  1974  -  January  1975, 
compares  the  results  of  the  wave  spectral  model  to  data  from 
the  NOAft.  Data  Buoy  Office  (NDBO)  EB-03  wave  bu'>y.  The  fourth 
case  demons tratas  an  application  of  the  wave  spectral  model. 

A  Ship  Routing  group  compared  output  of  the  spectral  model  to 
wave  observations  from  a  large  (950  feet)  container  ship.  The 


15 


container  ship  logged  wave  observations  and  these  observations 
were  compared  to  the  SOWM  analysis. 

A.  Case  Study  for  29  November-2  December  1969 

An  extratropical  storm  with  a  center  located  at 
approximately  40°N,  165°E  began  to  intensify  at  0600Z,  29 
November  1S69.  According  to  the  original  FNWC  wind  analysis 
the  maximum  wind  speed  attained  was  65  knots;  however,  the 
wind  fields  have  been  reanalysed  and  from  1200Z,  29  November 
60  0600,  1  December  the  highest  velocities  were  over  70  knots, 
with  a  maximum  wind  speed  of  91  knots  reached  at  1800Z,  30 
November.  The  storm  moved  in  a  southeasterly  direction  and 
the  wind  directions  associated  with  the  highest  speeds  ranged 
from  250°  to  340°.  The  wind  field  analysis  for  ]8Z,  30  November 
is  shown  in  Figure  13.  Only  every  other  grid  point  is  plotted; 
therefore,  the  highest  speed  may  not  be  shown.  The  highest 
wind  velocities,  their  location,  and  associated  significant 
wave  heights  (the  average  height  of  the  highest  one- third 
waves — see  Appendxx  A  Equation  12)  computed  by  the  SOWM  for 
every  six  hours  are  shown  in  Table  3. 

Using  the  original  wind  data  as  input,  the  FNWC  singular 
wave  model,  which  was  the  operational  wave  model  in  1969,  com¬ 
puted  a  maximum  significant  wind  wave  height  of  40  feet  (Hubert 
and  Mendenhall,  1970) .  This  wave  height  is  much  less  than  the 
significant  wave  heights  computed  by  the  wave  spectral  model 
(Table  3) .  Even  if  a  wind  speed  of  91  knots  had  been  used  as 
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input  to  the  singular  wave  model,  the  highest  significant  wind 
wave  height  computed  by  the  model  still  would  not  have  exceeded 
41  feet.  These  figures  demonstrate  one  of  the  limitations  of 
the  FNWC  singular  model. 

During  this  time  period  the  manned  spar  buoy  FLIP  was 
operating  in  a  vertical  mode  in  the  vicinity  of  27°30,N, 
157°45'W.  The  singular  wave  model  indicated  that  the  swell 
or  low  frequency  wave  energy  which  propagated  from  the  above- 
mentioned  storm  area  to  the  area  in  which  FLIP  was  located  had 
an  average  height  of  12-15  feet.  The  maximum  significant 
height  computed  by  the  singular  model  for  any  grid  point  during 
this  time  period  never  exceeded  25  feet.  Because  FLIP  was 
operating  in  a  vertical  mode,  scientists  aboard  FLIP  were 
able  to  accurately  measure  the  wave  heights  against  various 
points  on  the  buoy.  In  addition,  16mm  movies  were  taken  of 
FLIP  from  it's  tew  ship  and  the  data  from  the  movies  were 
correlated  with  the  FLIP  observations  (Rudnick  and  Hasse,  1971) . 
FLIP 1 s  log  indicated  that  the  seas  were  constantly  building  for 
approximately  35  hours  preceding  2300Z,  1  December.  The  period 
of  maximum  wave  intensity  existed  from  approximately  2300Z, 

1  December  to  0400Z,  2  December.  During  this  time  period  the 
wave  measurements  shown  in  Table  4  were  obtained.  It  is  obvious 
that  when  these  data  are  compared  to  the  singular  wave  model 
calculations,  the  singular  model  did  not  properly  propagate  the 
swell . 
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The  SOWM  wave  spectral  model  was  run  in  the  hindcast  mode 
for  the  time  period  of  16  November  -  3  December  1969.  The  model 
was  initialized  on  16  November  to  ensure  that  steady  state 
conditions  would  be  reached  by  29  November.  Plots  of  one¬ 
dimensional  wave  spectral  (spectral  density  vs.  frequency)  for 
the  closest  grid  point,  27.16°N,  156. 6°W,  to  FLIP'S  position 
are  shown  in  Figure  14.  The  wave  spectra  were  plotted  every 
six  hours  for  the  time  period  of  12Z,  01  December  to  06 Z,  2 
December.  Analyzed  wind,  velocities  are  also  indicated.  v The’ 
logs  of  near-by  ships  indicated  that  the  winds  for  this  time 
period  were  generally  from  the  west  and  had  speeds  on  the  order 
of  20-30  knots  (Rudnick,,et  al) .  Other  ship  logs  from  this 
area  (McConathy,  1974)  do  show  that  the  wind  speeds  started 
to  decrease  no  later  than  0930Z,  2  December.  Thus,  there 
appears  to  be  good  agreement  between  the  reported  and  analyzed 
wind  velocities  in  the  vicinity  of  FLIP. 

The  computer  wave  spectra  for  18Z,  30  November  and  0Z, 

6Z,  1  December  consist  primarily  of  high  frequency  wave  energy 
which  was  a  function  of  the  local  wind.  The  12Z  and  18  Z, 

1  December  wave  spectral  contained  both  high  and  low  frequency 
wave  energy.  The  last  two  wave  spectra,  0Z  and  6Z,  2  December, 
of  this  time  sequence  were  dominated  by  low  frequency  wave 
energy  whose  significant  frequency  bands  were  centered  at  0.04411" 
(22.7  seconds)  and  0,05HZ  (20  seconds).  As  shown  in  Table  4, 
the  observed  significant  frequency  bands  as  computed  from  the  movies 
which  were  taken  approximately  23Z  1  December  1969,  range  from  0.05H; 


to  0.062  HZ  (20  to  16  seconds).  The  entire  film  was  shot  during 
a  five  minute  interval  and  divided  into  five  sequencies  of  8, 

57/  45,  20  and  28  seconds.  Only  in  three  of  these  sequences 
would  there  have  been  enough  time  to  observe  wave  frequencies  of 
less  than  0.05  seconds  and  even  these  time  sequences  were  rather 
short.  It  is  not  unreasonable  to  assume  that  sometime  during 
the  five  hour  period  from  23Z,  1  December  to  04Z,  2  December, 
significant  wave  frequencies  of  less  than  0.05  HZ  could  have 
existed. 

Several  discrepancies  were  revealed  when  comparing  the 
analyzed  wave  data  to  the  measured  wave  data.  The  computed 
wave  energy  shown  in  Figure  14  peaked  at  06Z,  2  December  which 
was  at  least  two  hours  after  the  observed  wave  conditions 
actually  reached  their  greatest  intensity.  In  addition,  the 
computer  model  calculated  a  maximum  significant  wave  height 
of  49  feet.  Several  questions  arise:  (a)  are  these  discrep¬ 
ancies  solely  a  function  of  the  wave  spectral  model  or  are  there 
other  reasons  for  the  errors. 

In  this  case  the  SOWM  analysis  were  far  superior  to  the 
singular  model  analysis.  Had  the  SOWM  been  operational  in  1969, 
(even  with  discrepancies  described  above)  FLIP  could  have  been 
better  prepared  for  the  intense  low  frequency  wave  energy  which 
propagated  through  it's  operating  area.  It  is  interesting  to 
note  that  FLIP  was  able  to  withstand  wave  heights  up  to  80  feet 
with  a  minimum  of  damage;  however,  as  shown  in  Figure  15  (Kerr, 
1964),  FLIP  was  probably  in  more  danger  because  the  significant 
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frequency  of  the  propagating  wave  energy  was  approaching  the 
resonant  frequency  of  FLIP.  If  the  significant  frequency  of 
the  wave  energy  had  equaled  the  resonance  frequency  while  FLIP 
was  in  the  vertical  position#  FLIP  would  become  unstable  and 
probably  capsized.  The  SOWM  did  accurately  depict  the  range 
of  the  significant  frequency  bands  of  the  low  frequency  wave 
energy. 

The  wave  spectral  model  propagated  a  larger  portion  of 
the  low  frequency  wave  energy  to  the  north  of  the  FLIP.  This 
may  have  been  caused  by  errors  in  the  analyzed  wind  fields 
used  as  input  to  the  model.  If  the  storm's  geographical  po¬ 
sition  was  not  analyzed  correctly  or  the  wind  directions  with¬ 
in  the  storm  area  were  incorrect,  the  low  frequency  energy 
would  not  propagate  in  the  right  direction.  The  low  frequency 
wave  energy  in  the  area  of  FLIP  seemed  to  be  concentrated  in 
a  rather  narrow  direction  band  as  can  be  seen  in  Figures  16  and 
17.  Prior  to  C6Z,  2  December  almost  all  the  low  frequency  wave 
energy  was  concentrated  in  the  308°  direction  band  and  at  06Z# 

2  December  the  wave  energy  was  primarily  located  in  the  308° 
and  338°  direction  bands.  The  grid  points  to  the  north  of 
FLIP  also  had  significant  wave  energy  in  the  278°  and  248° 
bands  as  seen  in  Figure  18.  The  significant  wave  height  for 
the  grid  point  (39 °N,  161°W)  shown  in  this  figure  was  com¬ 
parable  to  that  measured  by  FLIP.  Since  there  was  no  measured 
data  available  from  the  storm  area,  it  is  possible  that  the 
analysis  incorrectly  positioned  the  storm  or  miscalculated  the 
wind  direction. 
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This  case  study  does  point  out  the  type  of  errors  that  can 
occur  in  the  wave  spectral  model  if  the  wind  fields  are  in¬ 
correct.  An  error  of  15“  in  the  wind  direction  can  be  signifi¬ 
cant  to  the  wave  spectral  model  since  the  wave  energy  is  divided 
into  30°  direction  bands.  If  the  wave  energy  is  misplaced  by 
one  direction  band,  and  this  wave  energy  is  allowed  to  propagate 
with  minimum  interference,  the  ensuing  can  cause  the  wave  energy 
to  be  misplaced  by  several  grid  points  over  a  1000  mile  track. 
This  could  have  occurred  in  this  case  study. 

In  addition  to  improving  the  wind  analysis,  the  above  prob¬ 
lem  might  have  been  partially  alleviated  by  increasing  the  number 
of  direction  bands  used  to  define  the  wave  energy  spectra.  Pro¬ 
fessor  Pierson  divides  the  spectra  into  15®  increments  for 
hindcast  studies.  Of  course,  the  number  of  parameters  to  be 
calculated  at  each  grid  point  increases  from  180  to  360,  re¬ 
quiring  additional  computer  memory  and  computation  time  to  run 
the  model,  currently  impractical  on  the  FNWC  computer  system. 

B.  Case  Study  -  26  October  -  28  October  1973 

Until  December  1974  wave  spectra  computed  from  wave 
recorders  on  a  routine,  synoptic  time  basis  were  not  available. 
Wave  recorders  usually  were  installed  for  specific  projects 
and  after  the  projects  were  completed,  the  recorders  were  re¬ 
moved  from  their  moorings.  To  the  best  of  our  knowledge  there 

I 

has  only  been  one  deep-water  wave  recorder  being  operated  in  a 
semi-permanent  basis  in  the  North  Pacific  Ocean.  The  Canadians 
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keep  a  weather  ship  at  50°N,  145°W,  Ocean  Station  Vessel  (OSV) 
PAPA,  as  seen  in  Figure  19  and  the  ships  which  maintain  this 
station  have  wave  recorders  mounted  in  their  bow. 

On  26  and  28  October  1973  the  Canadian  weather  ship,  CCG&S 
VANCOUVER  ,  made  special  wave  measurements  in  conjunction 
with  a  NAVOCEANO  wave  project  at  OSV  PAPA.  Both  sets  of 
measurements  began  at  approximately  20Z  and  lasted  for  about 
one  hour.  These  measurements  were  made  in  addition  to  the  nor¬ 
mal  wave  measurements  madei  every  three  hours  at  OSV  PAPA  and 
were  probably  of  better  quality  than  the  routine  measurements 
since  the  ship  was  oriented  in  the  best  position  to  record  the 
waves.  - 

The^ Tucke.r  shipboard ’wave  recorders  used  by  the  Canadian  “ 
weather  ships  have  been  described  by  Tucker  (1956)  and  only 
significant  points  will  be  discussed  in  this  paper.  The  Pierson' 
Moskowitz  fully  developed  spectra  were  derived  from  Tucker  wave 
data.  The  pressure  transducer  measures  the  high  frequency  wave 
energy  up  to  0.15HZ.  This  upper  limitation  is  due  to  the 
attentuation  of  the  wave  pressure  with  depth  (Briscoe,  1971) . 

The  wave  recorders  are  mounted  approximately  13  feet  below  the 
mean  water-line  on  the  Canadian  weather  ships.  The  vertical 
accelerometer  measures  the  low  frequency  energy  down  to  0.04HZ 
and  at  that  frequency,  errors  in  the  electronic  double  inte¬ 
gration  of  the  accelerometer  signal  used  to  cancel  ship's  heave 
becomes  significant.  Since  the  frequency  boundaries  are  essen¬ 
tially  0.16KZ  and  0.039HZ  in  the  wave  spectral  model,  the 
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frequency  range  of  the  Tucker  recorder  is  quite  adequate.  The 
ship  must  be  in  a  hove-to  position,  preferably  headed  into  the 
waves,  in  order  to  accurately  measure  the  wave  heights.  The 
data  from  the  Tucker  recorder  are  currently  recorded  as  analog 
traces  by  a  strip  chart  recorder  at  a  speed  of  either  one  inch 
or  2  inches  per  minute  for  approximately  15  to  20  minutes  every 
three  hours.  The  analog  traces  were  digitized  every  0.01  inches 
which  is  equivalent  to  0.6  or  0.3  seconds,  depending  on  the 
analog  chart  speed,  at  FNWC.  Before  wave  spectral  could  be 
computed  from  the  digitized  wave  data,  a  correction  factor  was 
applied  to  the  wave  data.  Several  versions  of  the  correction 
factor  have  been  cited  in  the  literature  (M.  Darbyshire,  1960; 

J.  Darbyshire,  1963;  D.  Cartwright,  1963  Ewing,  1969;  Holland, 
1972);  however,  the  most  common  correction  is  shown  to  be: 

H 

JL  .  03  {  w  )  A  (K)  (1) 

Hr 

and  C  is  a  constant  (0.83) 

B  (  iti )  is  due  to  the  characteristics  of  the  internal  electri¬ 
cal  circuits  of  the  instrument  and  is  equal  to  [1  +  (8.8  w)“2]3/‘ 

where  w  »  2  r/T  ,  T  is  the  near  crossing  period  (i.e.,  the 
in  m 

number  of  positive  zero  crossing  divided  by  the  total  time  and 
8.8  represents  the  time  constant  of  three  filters  used  in  the 
instrument. 

A  (K)  represents  the  relation  between  the  height  of  water 
above  a  fixed  point  on  the  ship's  hull  and  the  pressure  measured 


at  this  point.  A(K)  is  also  affected  by  the  ship's  heaving 
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^  response  and  is  equal  to  e2*5du)  where  d  is  the  depth  of 
the  recorder  below  the  mean  water  line  (13  feet)  and  g  is  the 
acceleration  of  gravity.  The  ratio  of  the  true  wave  height 
(Ht)  to  the  measure  wave  height  (Hr)  can  be  written: 


B*  =  0.83  (1  +  (8.8  *r2i3/2  e-2'5  (2» 

r 

Note  that  the  correction  factor  is  a  function,  of  frequency  and 
the  higher  the  frequency  the  larger  the  correction  factor. 

Table  5  shows  the  ratios  that  can  be  applied  depending  on  the 
mean  cross  frequency  (period).  Dr.  Briscoe  (1971)  suggested 
that  no  correction  factor  be  applied  to  the  wave  records.  He 
arrived  at  this  decision  based  on  information  supplied  to  him  by 
the  National  Institute  of  Oceanography.  Briscoe  believes  that 
it  is  difficult  to  select  one  correction  factor  which  is  depen¬ 
dent  on  all  frequencies  in  the  spectrum  and  there  is  less  overall 
error  whar.  no  correction  factor  is  applied  than  when  an  incorrect 
factor  is  used. 

The  Canadians  alternate  two  weather  ships  at  OSV  PAPA.  Dr. 
J.  R.  Wilson  of  the  Canadian  Marine  Environmental  Data  Service, 
who  graciously  allowed  us  to  use  the  original  analog  wave 
records,  did  some  comparison  between  corrected  data  from  the 
Tucker  recorders  used  aboard  the  CCG  t  S  QUADRA  and  from  a  Data- 
well  Wave  rider  (personal  communication,  1974) .  The  two  sets 
of  data  compared  rather  well.  Data  from  the  CCG  &  S  VANCOUVER 
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have  not  been  compared  to  data  from  other  instruments;  however, 
there  is  a  possibility,  based  on  hand  analyses,  that  the  cor¬ 
rected  wave  heights  from  this  recorder  may  be  too  high. 

Table  6  contains  the  observed  wind  velocities  from  OSV 
PAPA  and  the  analyzed  wind  velocities  from  the  two  closest  grid 
points  to  OSV  PAPA  for  the  time  period  of  06Z,  26  October  -  OZ, 

29  October.  The  analyzed  wind  fields  were  computed  every  six 
hours  while  the  observed  winds  were  obtained  every  three  hours. 
Generally  speaking,  the  analyzed  winds  compared  rather  well 
with  the  measured  winds.  In  some  cases  such  as  18Z,  26  October 
the  wind  velocities  for  one  grid  point  may  compare  better  with 
the  observed  wind  velocities  than  the  wind  velocities  from  the 
adjacent  grid  point.  These  differences  will  be  reflected  in 
the  computation  of  the  wave  spectra. 

At  06Z,  24  October  there  was  a  low  pressure  center  located 
at  50°N,  178W,  approximately  1000  nautical  miles  from  OSV  PAPA. 
Initially,  the  highest  wind  speeds  around  the  low  pressure 
center  were  on  the  order  of  30-35  knots,  but  as  the  low  pressure 
center  moved  in  a  northeasterly  direction,  the  pressure  gradients 
became  less  intense  and  the  highest  wind  speeds  decreased  to 
20-25  knots.  A  plot  of  the  12Z,  25  October  wind  field  is  shown 
in  Figure  19.  It  can  be  seen  from  this  figure  that  the  wave 
energy  propagated  from  this  low  pressure  area  towards  OSV  PAPA 
came  from  approximately  200°-250°.  At  12Z,  26  October  the 
analyzed  low  pressure  center  was  to  the  northwest  of  OSV  PAPA 
and  the  pressure  gradients  around  the  center  were  intensifying. 
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A. 


The  analyzed  wind  field  for  122,  26  October  is  shown  in  Figure  20 

Comparisons  between  non  directional  wave  spectral  computed 

from  OSV  PAPA  data  and  the  2  grid  points  of  the  SOWM  are  shown 

in  Figures  21-28  for  18Z,  20Z,  26  October  and  002,  27  October. 

The  wave  spectral  model  computes  wave  spectra  in  terms  of  var- 
2 

iance  (ft  ) ;  however,  since  the  fifteen  band  widths  are  unequal, 
it  was  decided  the  wave  spectra  would  be  better  represented  by 
variance  density  {ft  /HZ)  which  is  obtained  by  multiplying  the 
variance  by  the  frequency  band  widths.  There  are  seme' instances 
such  as  an  example  which  will  be  shown  later  where  this  process 
may  obscure  spectral  peaks.  Even  when  the  wave  spectra  are 
converted  to  variance  density,  there  are  instances  when  the 
curves  -are  not  completely  smoothed.  For  example, '  the" double 
peaks  in  Figures  22,  24,  25  and  28  are  probably  a  function  of 
the  frequency  band  widths  and  do  not  really  exist. 

The  OSV  PAPA  wave  spectra  were  computed  by  the  Blackman- 
Tukey  method  (1959’ .  Thirty  degrees  of  freedom  was  selected  for 
the  computation  cf  all  the  spectra.  This  means  that  the  number 
of  lags  varied  for  each  spectrum,  depending  on  the  length  of 
wave  record  that  was  digitized. 

Generally  speaking,  the  computed  wave  spectra,  including 
significant  wave  heights  and  dominant  frequencies,  compared 
favorably  with  the  measured  wave  spectra.  In  some  case  such 
as  shown  in  Figures  24  and  25  the  comparison  was  better  for 
one  grid  point  (in  this  case,  48.5°N,  142. 5°W)  than  for  the 
adjacent  grid  point  (50.9*N,  145. 6°W).  In  this  instance  the 
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discrepancy  can  be  accounted  for  by  the  differences  in  the 
analyzed  and  observed  wind  speeds.  The  significance  of  accurate 
wind  velocity  input  can  be  seen  in  Figure  25  where  the  computed 
wave  spectrum  for  OZ,  26  October  is  compared  to  the  measured 
wave  spectrum  for  20Z,  26  October  and  Figure  27  where  the  same 
computed  wave  spectrum  is  compared  to  the  measured  spectrum  for 
OZ,  27  October.  The  comparison  is  better  than  in  Figure  27 
because  a  higher  wind  speed  existed  in  OSV  PAPA  for  at  least 
three  more  hours  than  was  analyzed  at  the  grid  point.  The  low 
frequency  energy  of  the  computer  wave  spectra  for  48.5°N, 

140. 5°W  compares  better  to  the  low  frequency  energy  that  was 
measured  at  OSV  PAPA  (Figures  24  and  28) .  This  also  can  be 
seen  in  Figures  29  and  30  which  contain  the  computed  two  di¬ 
mensional  wave  spectra.  More  low  frequency  energy  exists  at 
48 . 5°N,  152. 5°W  than  at  50.9°N,  145. 6°W  which  is  closer  to 
OSV  PAPA.  Note  that  at  18Z,  26  October  (Figure  29)  the  pre¬ 
dominant  low  frequency  energy  has  shifted  to  200°  at  both  grid 
points.  Thus  it  appears  that,  in  this  case,  given  correct 
analyzed  winds,  the  computer  model  is  propagating  the  low  fre¬ 
quency  energy  correctly. 

A  rather  poor  comparison  at  06 Z,  27  October  is  shown  for 
48 . 5°N,  142. 5°W  in  figure  31.  The  comparison  for  the  adjacent 
grid  point  is  similar.  The  computed  spectrum  compares  favorably 
to  the  measured  spectrum  in  the  low  frequency  range  but  does  not 
compare  at  all  in  the  high  frequency  range.  An  investigation  of 
the  computed  two  dimensional  wave  spectrum  which  is  calculated 
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in  terms  of  variance  indicates  that  an  energy  peak  existed  at 
0.092HZ.  The  energy  peak  was  obscured  when  the  spectrum  was 
converted  to  variance  density  by  multiplying  by  the  uneven 
bandwidths.  The  measured  energy  peak  which  was  significantly 
larger  in  magnitude  was  at  0.105HZ.  A  similar  situation  occurred 
between  18Z,  28  October  and  OZ,  29  October  and  will  be  discussed 
in  a  later  section.  Some  possible  causes  of  these  discrepancies 
will  also  be  discussed. 

This  case  indicates  .the  importance  of  being  aware  of  which 
physical  units,  i.e.,  variance  or  variance  density,  are  being 
used  to  display  the  model  wave  spectra.  Since  the  frequency 
bandwidths  are  unequal,  two  differently  shaped  curves  can.be 
obtained  from  the  same  set  of  data.  For  example,"  if  only  the 
variance  density  wave  spectra  were  available  for  this  case,  the 
high  frequency  energy  peak  would  have  been  obscured  by  the  low 
frequency  energy;  however,  a  review  of  the  variance  wave  spectra 
indicates  that  there  was  significant  high  frequency  wave  energy 
at  the  grid  point. 

After  00Z,  27  October  the  low  pressure  center  mentioned 
earlier  continued  to  move  northeastward.  The  12Z,  27  October 
analysis  places  the  low  pressure  center  directly  to  the  north 
of  OSV  PAPA  (Figure  32).  At  18Z,  27  October  a  new  low  pressure 
center  formed  southwest  of  OSV  PAPA  (Figure  33) .  The  wind 
system  which  generated  the  wave  energy  that  propagated  towards 
OSV  PAPA  had  speeds  of  30-40  knots  from  the  south  (Figure  34) . 
Thus,  between  18Z,  28  October  and  OZ,  29.  October  two  distinct 
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wave  trains  existed  at  .OSV  PAPA  which  were  confirmed  by  photo¬ 
graphs  taken  by  a  NAVOCEANO  airplane  between  20-21Z,  28  October. 

There  were  no  wave  spectra  from  OSV  PAPA  available  for 
comparison  between  00Z,  28  October  and  2QZ,  28  October.  The 
18Z,  28  October  and  00Z,  29  October  computer  wave  model  wave 
spectra  for  the  grid  point  at  5Q.9°N,  145. 6°W  are  compared 
to  the  20Z  and  21Z  wave  spectra  from  OSV  PAPA  in  Figures  35- 
37.  The  OSV  PAPA  wave  spectra  were  computed  from  twenty  minute 
records  for  each  hour.  The  OSV  PAPA  spectra  in  Figures  35-36 
indicate  that  the  wave  energy  was  growing  rapidly  as  a  function 
of  the  local  wind  velocity.  The  high  frequency  wave  energy 
shifted  towards  the  low  frequency  end  of  the  scale.  Data  ob¬ 
tained  by  a  NAVOCEANO  airplane  during  the  same  time  period 
agreed  with  the  CSV  PAPA  spectra.  The  computer  model  wave 
spectrum  also  had  two  peaks  expected,  but  as  in  the  OGZ,  27 
October  case,  the  high  frequency  energy  peak  is  located  at  a 
lower  frequency  than  the  OSV  PAPA  high  frequency  energy  peaks 
even  though  the  CSV  PAPA  wave  spectrum  in  Figure  36  was  measured 
three  hours  after  the  model  wave  spectrum  was  computed.  Un¬ 
fortunately,  the  weather  ship  was  underway  at  OZ,  29  October  and 
the  measured  wave  spectrum  obtained  at  that  time  was  considered 
unreliable.  The  analyzed  wave  spectrum  as  shown  in  Figure  37 
no  longer  has  bi-nodal  peaks  and  the  wave  energy  is  concentrated 
in  the  lower  frequency  range. 

A  review  of  the  analyzed  and  observed  wind  velocities  listed 
in  Table  6  for  both  06Z,  27  October  and  18Z  October-OZ,  29  October 
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seems  to  indicate  that  the  two  sets  of  winds  were  in  good  agree¬ 
ment.  In  the  majority  of  the  comparisons  the  analyzed  wind  speeds 
were  lower  than  the  observed  wind  speeds.  The  one  significant 
difference  between  the  observed  and  analyzed  wind  velocities 
is  that  the  observed  winds  show  change  every  three  hours  while 
the  analyzed  winds  which  are  computed  every  six  hours,  or  two 
wave  model  time  steps.  Thus,  the  12z,  28  October  wind  speed  of 
21.8  knots  was  used  as  input  at  09Z  and  12Z  and  the  18Z,  28 
October  wind  speed  of  34.3  knots  was  used  at  15Z  and  18Z.  The 
differences  between  the  observed  and  analyzed  wind  speeds  at 
09Z  and  15Z  were  approximately  8  knots  for  each  case.  This 
means  that  there  were  two  time  steps  over  a  twelve  hour  period 
where  the  computed  wave  energy  was  being  grown  as  a  function  of 
a  significantly  higher  wind  velocity  than  what  actually  existed 
at  that  time;  however,  this  argument  can  not  be  used  for  the 
06Z,  27  October  case  wnen  the  analyzed  wind  speeds  at  50.9°N, 

145. 6°W  were  always  less  than  the  observed  wind  speeds. 

The  wind  directions  in  both  cases  were  such  that  if  the 
energy  dissipation  equation  (Appendix  A,  equation  11)  in  the 
computer  model  is  correct,  energy  dissipation  due  to  wind 
velocity  should  be  a  minimum  during  this  time  period.  An 
investigation  of  the  wind  directions  do  point  out  one  possible 
explanation  as  to  what  may  have  caused  the  discrepancies  in  the 
high  frequency  range.  In  both  the  06Z,  27  October  and  18Z, 

28  October-OCZ ,  29  October  cases,  the  wind  was  primarily  from 
the  south  to  southwest  for  at  least  six  hours  prior  to  and  at 
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analysis  time.  A  southerly  wind  should  increase  the  air  temper¬ 
ature.  If  the  air  temperature  was  significantly  higher  than  the 
underlying  sea  surface  temperature,  then  stable  atmospheric 
conditions  would  exist.  Cardone  has  used  the  air-sea  temperature 
difference  of  0.9°C  as  the  cut-off  point.  As  mentioned  in 
Section  III,  the  wind  velocities,  which  are  derived  from  the  wind 
fields  and  used  as  input  to  the  growth  portion  of  the  wave  spectral 
model,  are  calculated  for  neutral  stability.  If  stable  conditions 
existed,  the  wave  energy  growth  rate  in  the  model  may  have  been 
excessive.  A  review  of  ship  observations  for  the  area  during 
both  time  periods  could  neither  confirm  or  disprove  this  hypoth¬ 
esis  as  temperature  data  from  several  ships  in  the  area  varied 
considerable. 

Another  possibility  is  the  growth  equations  and/or  the 
Pierson-Moskowitz  (P/M)  wave  spectra  used  in  the  computer  model 
do  not  always  represent  the  existing  sea  conditions.  A  similar 
hypothesis  has  been  suggested  by  W.  E.  Cummins  (1974)  of  the 
Naval  Ship  Research  and  Development  Center.  Cummins  compared 
P/M  response  curves  to  wave  spectral  from  OSV  INDIA  (60°N,  20°W) 
while  investigating  ship  responses  in  various  sea  conditions. 

OSV  PAPA  and  OSV  INDIA  have  several  points  in  common.  Since 
extratropical  storms  propagate  primarily  from  west  to  east, 

OSV  INDIA  which  is  located  in  the  northeastern  Atlantic  is  able 
to  measure  wave  energy  across  a  broad  frequency  range,  as  does 
OSV  PAPA  ir.  the  northeastern  Pacific  Ocean. 

Extratropical  storms  in  the  Pacific  Ocean,  generally  speaking, 
are  probably  more  intense  and  cover  a  larger  area  than  storms  in  the 
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Atlantic  and  wave  spectra  in  the  Pacific  probably  cover  a  larger 
frequency  range  and  contain  more  wave  energy.  Both  OSV  PAPA  and 
OSV  INDIA  use  Tucker  shipboard  recorders  to  measure  waves  and 
the  P/M  curves  were  derived  from  Tucker  wave  recorder  data; 
thus,  all  of  these  spectra  should  have  the  same  basic  charac¬ 
teristics.  Cummins  believed  that  since  P/M  spectra  essentially 
represent  a  theoretical  sea;  i.e.,  generated  by  a  steady  wind 
blowing  for  an  infinite  time  over  an  infinite  fetch,  uncontam- 
minated  by  swell,  these  spectral  will  fall  to  one  side  of  the 
OSV  INDIA  spectra  because  the  measured  spectra  reflect  the 
total  energy — variable  locally  generated  wave  energy  and  energy 
propagated  from  distant  generating  areas — at  the  weather  station. 
Cummins  cited  three  sea  conditions  where  the  P/M  curves  had 
major  discrepancies:  in  low  significant  wave  height  situations 
where  the  swell  dominated  the  P/M  spectra  were  too  low  ;in  the 
intermediate  wave  heights  cases  the  P/M  spectra  were  significantly 
larger  than  the  OSV  INDIA  spectra  which,  Cummins  claimed,  meant 
that  these  seas  were  rarely  fully  developed;  in  the  highest  in 
wave  height  conditions  the  P/M  spectra  were  low,  sometimes  by 
a  factor  of  two,  which  Cummins  suggest  meant  these  very  high 
seas  had  much  shorter  wave  lengths  than  the  P/M  spectra  would 
imply.  This  could  arise  from  very  strong  winds  blowing  for 
short  periods  (not  defined),  with  the  result  that  high  frequency 
steep  seas  would  be  generated  but  would  never  become  fully 
developed. 
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Since/  as  mentioned  in  Section  IIA,  the  icosahedral  gnomonic 
wave  spectral  model  uses  the  P/M  spectra  as  upper  limits  to  its' 
growth  curves  and  only  in  the  fully  developed  state  would  the 
model  spectra  resemble  the  P/M  spectra,  the  first  two  situations 
described  by  Cummins  may  not  occur;  however,  the  last  situation 
is  similar  to  the  two  cases  cited  above.  In  those  cases  there 
seem  to  be  two  distinct  wave  trains* — one  consisting  of  low  fre¬ 
quency  wave  energy  propagated  from  adjacent  points  and  the  other 
containing  energy  locally  generated.  This  can  be  seen  in  the 
OSV  PAPA  wave  spectra  (Figures  31,  35-37)  which  have  two  distinct 
energy  peaks.  The  model  spectra  only  have  low  frequency  peaks. 

If  Cummins'  conclusions  are  correct,  then  it  is  possible  that 
under  certain  wind  conditions  the  P/M  curves  are  restricting  the 
growth  rate  of  the  wave  energy  in  the  high  frequency  range  and 
forcing  the  wave  energy  to  shift  toward  the  low  frequency  range 
sooner  than  it  should.  Certainly  in  the  two  cases  under 

discussion  this  explanation  seems  plausible. 

More  data  are  needed  for  investigation  of  the  two  above 
hypotheses.  It  is  hoped  that  data  from  the  NDBO  wave  buoys 
will  provide  some  of  the  answers. 

C.  Case  Study  -  December  1974  -  January  1975 

A  NOAA  Data  Buoy  Office  (NDBO)  wave  buoy  (EB03)  was 
deployed  in  the  northeastern  Pacific  Ocean  (Gulf  of  Alaska) 
at  56.0°N,  147, 9°W  in  3  December  1974  in  a  location  where 
severe  extratropical  storms  transit  during  the  winter  season. 

Since  most  of  the  storms  develop  in  the  northwest  Pacific  Ocean, 
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this  buoy  is  ideally  located  to  monitor  the  propagation  of  wave 
energy  over  great  distances  as  well  as  every  growth  for  local 
high  wind  velocity  cases.  The  buoy  is  configured  to  provide 
wave  spectral  and  wind  information  in  real  time  to  a  shore 
station. 

The  buoy  itself  is  a  large  discus  (40  feet  in  diameter) 
which  NDBO  claims  exhibits  a  near  unity  response  amplitude 
ratio  over  a  broad  frequency  band  below  a  high  frequency 
cutoff  (Michelena,  Steele,  Niedermann  and  Hindman,  1974) .  ' 
Assuming  this  to  be  true, ^monitoring  the  buoy  motion  will  give 
a  response  equivalent  to  the  sea  surface  motion.  The  wave 
sensing  instrument  is  a.  servo- type  accelerometer  whose  axis 
is  fixed  to  the  vertical  center  line  of  the  buoy.  Since  the 
gyro-stabilized  platforms  which  NDBO  had  used  previously  to 
maintain  the  axis  of  the  linear  accelerometer  vertical  relative 
to  the  mean  surface  have  proved  to  be  unreliable,  NDBO  realized 
that  by  fixing  the  axis  of  the  accelerometer,  an  error  would  be 
introduced  and  that  the  determination  of  the  error  magnitude 
would  be  diffcult  because  of  the  interactions  of  the  buoy  hull 
and  the  sea  surface.  NDBO  made  the  assumptions  that  the  sea 
was  an  unidirectional  simple  swell  wave  of  one  frequency  and 
the  buoy  was  a  perfect  surface-follower  that  moves  as  a  particle 
would — in  a  circle.  Based  on  highly  idealized  sea  conditions 
and  buoy  motion,  NDBO  found  the  error  in  wave  heights  to  be  on 
the  order  of  24.  An  assumption  was  made  that  the  error  would 
not  be  significantly  larger  under  more  complex  conditions. 
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The  analog  wave  data  are  digitized  every  0.5  seconds  over  a 
fifteen  minute  period  (1800  points)  every  three  hours.  Wave 
spectra  are  computed  from  the  acceleration  data  and  are  con¬ 
verted  to  displacement  spectral  densities  by  multiplying  by 

4 

u>  (circular  frequency) .  The  spectra  are  considered  accurate 
within  a  frequency  range  of  0.01HZ  to  0.5HZ.  Significant  wave 
height,  average  wave  period  and  spectral  width  are  obtained  from 
moments  of  the  displacement  spectrum. 

The  wave  buoy  was.  calibrated  in  the  laboratory  comparisons 
were  made  between  data  obtained  at  sea  by  the  wave  buoy  and 
from  nearby  a  Datawell  Wave  Buoy  a  highly  reliable  and  accurate 
wave  sensor.  The  comparison  between  the  wave  spectra  from  the 
wave  sensor's  was  excellent.  A  report  has  been  written  by  Ocean 
Data  Systems,  Inc.  (ODSI)  for  NDBO  Comparing  wave  spectral  from 
the  wave  buoy  to  wave  spectra  from  the  SOWM  (1975) .  One  of  the 
conclusions  of  this  report  was  that  the  significant  wave  heights 
computed  from  the  SOWM  were  generally  higher  than  the  significant 
wave  heights  computed  from  the  buoy  data.  ODSI  concluded  that 
the  computer  model  wave  spectra  had  20%  excess  energy  and 
suggested  that  the  cause  was  lack  of  strong  decay  coefficients 
in  the  low  frequency  range.  Another  conclusion  was  that  the 
SOWM  in  some  cases  seemed  to  anticipate  the  arrival  of  the 
low  frequency  wave  energy.  Although  significant  if  true, 
these  conclusions  were  apparently  made  by  only  comparing 
significant  wave  heights,  average  wave  periods,  and  wave  spectra, 
and  without  regard  to  the  meteorological  conditions  either  up¬ 
stream  or  at  the  grid  point.  Without  further  analysis,  it 
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is  difficult  to  ascertain  whether  the  errors  are  a  function  of 
the  wind  input  or  the  wave  spectral  model  itself.  Some  of  the 
comparative  spectra,  along  with  the  local  analyzed  and  observed 
wind  velocities,  are  shown  in  Figures  38-42.  Most  of  the  com¬ 
parisons  are  good  and  many  of  the  differences  seem  to  be  caused 
by  differences  in  the  local  wind  velocities.  In  one  case  (Figure 
42)  where  there  was  low  velocity  winds  there  was  significant  wave 
energy  fn  the  frequency  bands  centered  at  approximately  0.08HZ 
of  both  spectra  that  was  either  propagated  from  other  areas  or 
had  been  generated  before  the  local  wind  had  subsided.  In 
several  of  the  wave  spectra  not  presented  in  this  paper  there 
does  seem  to  be  differences  which  can*t  be  accounted  for  by 
local  wind  velocity  differences.  Since  the  wave  spectral  data 
from  EB03  are  being  transmitted  routinely  to  FNWC,  continuing 
analysis  and  comparison  of  the  data  are  being  made. 

D.  Case  Scucy  -  Application  of  Wave  Spectral  Data  to 
Optimum  Track  Ship  Routing 

One  of  the  principal  applications  of  wave  spectral 
data  is  Optimum  Track  Ship  Routing  (OTSR) .  There  are  two 
U.S.  Navy  OTSR  groups — one  of  which  is  located  at  FNWC  and 
responsible  for  U.S.  Navy  ships  and  ships  under  military  con¬ 
tract  that  transit  the  Pacific  and  Indian  Oceans  and  the  other, 
located  at  FWC,  Norfolk  is  responsible  for  ships  that  transit 
the  Atlantic  Ocean.  It  has  been  demonstrated  that  the  OTSR 
groups  save  the  U.S.  Navy  at  least  4.5  million  dollars  per 
year  by  minimizing  fuel  consumption  during  transit,  minimizing 
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ship  damage  in  storm  areas,  and  increased  reliability  of 
ship ' s  schedules . 

A  container  ship's  log  containing  positions,  wave  period, 
height  and  direction  observations,  and  ship's  speed  made 
while  the  ship  transitted  the  Pacific  Ocean  has  been  obtained. 
Although  visual  observations  are  usually  considered  unreliable, 
observations  from  this  ship  have  proved  to  be  highly  reliable 
in  the  past  and  in  this  case  were  made  by  the  Master  in  reponse 
to  a  request  for  carefully  observed  data.  The  ship  is  950  feet 
in  length  and  capable  of  speeds  in  excess  of  30  knots.  Excerpts 
of  the  ship's  log  compared  to  output  of  the  computer  wave 
spectral  model  are  shown  in  Table  7.  The  actual  ship’s  track 
along  with  the  significant  wave  height  analysis  for  00Z,  7 
February  1975  are  shown  in  Figure  43.  Considering  that  the 
analyzed  wave  information  in  Table  7  were  only  available  on 
twelve  hour  increments,  the  analyzed  wave  heights  and  directions 
compared  very  well  to  the  observed  data.  The  wave  directions 
shown  in  parentheses  are  the  secondary  directions.  The  most 
interesting  comparisons  occurred  between  22Z,  6  February  and 
05Z,  7  February.  Prior  to  this  time  the  ship  had  maintained 
speeds  on  the  order  of  26-27  knots.  At  22Z,  6  February  the 
ship's  speed  began  to  decrease  and  for  several  hours  (there 
was  no  notation  in  the  log  from  05Z-21Z,  7  February)  the  ship 
was  unable  to  maintain  a  speed  greater  than  8  knots  because 
of  high  wave  conditions.  At  00Z,  7  February  the  ship's  ob¬ 
servations  was  50  feet  while  the  SOWM  had  calculated  a  wave 
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height  of  42  feet.  The  wave  directions  also  compared  well. 

*• 

It  can  be  seen  from  Figure  35  that  the  model  did  compute  sig¬ 
nificant  wave  heights  up  to  55  feet.  These  waves  were  approx¬ 
imately  300  nautical  miles  northwest  of  the  ship's  position. 

Based  on  the  distance  the  ship  travelled  from  05Z-21Z,  it  appears 
that  the  ship  could  not  average  more  than  14  knots  during  this 
time  period.  In  order  to  make  up  the  lost  time  the  ship  had 
to  maintain  speeds  of  30-32  knots  for  the  remainder  of  the 
journey  which  required  greater  fuel  consumption  and  thus,  greater 
operating  costs  for  the  transit. 

It  is  interesting  to  note  that,  based  on  the  ship's  log, 
the  ship  could  maintain  normal  operating  speeds  in  wave  heights 
up  to  25  feet. 

Prior  to  the  advent  of  the  SOWM,  OTSR  used  hand  analyzed 
wave  charts  and  the  output  of  the  singular  wave  model  was  used 
to  route  ships.  The  FNWC  singular  wave  model  had  a  limiter 
which  kept  it  from  building  wave  heights  greater  than  about 
44  feet  and  since  the  constraint  conditions  for  OTSR  routes 
and  diversions  were  wave  heights  of  20  feet  or  less  the  OTSR 
were  frequently  net  carefully  analized  in  the  areas  where 
sea  heights  were  extremely  high.  The  wave  product  users  have 
not  heretofore  beer,  used  to  seeing  analyzed  significant  wave 
heights  of  43  or  greater.  It  has  been  demonstrated  that  the 
SOWM  is  doing  a  better  job  of  representing  the  total  range  of 
energy  in  the  ocean  than  the  singular  model  or  hand  analysis 
procedures  did. 
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VI .  CONCLUSION 


The  operational  Icosahedral-Gnomonic  wave  spectral  com¬ 
puter  model  (SOWM) ,  produces  a  far  superior  wave  product  than  the 
FNWC  singular  wave  model/  the  previous  operational  model.  The 
singular  wave  model  was  limited  from  growing  waves  greater 
than  44  feet  and  could  not  adequately  describe  complex  wave 
conditions.  The  SOWM  considers  the  total  wave  energy  in  a  15 
frequency  band  by  12  direction  bands  matrix  in  a  grid  point 
and  propagates  the  energy  throughout  the  grid  system  as  a 
function  of  frequency  and  direction.  The  case  studies  in 
Section  V  demonstrated  the  ability  of  the  wave  model  to  accur¬ 
ately  propagate  low  frequency  energy  (Study  A) ,  to  cope  with 
more  than  one  wave  train  at  any  given  time  (Study  B) ,  to 
correlate  well  with  observations  from  an  operational  spectral 
observing  station  (Study  C) ,  and  to  grow  high  waves  for  high 
wind  conditions  fStudy  D) . 

Since  this  is  the  first  time  that  a  hemispheric  wave 
spectral  model  has  been  used  routinely  on  an  operational  basis, 
it  is  expected  that  some  modifications  will  have  to  be  made 
to  improve  the  model.  Other  wave  spectral  models  such  as  the 
Gelci  spectrcangular  model  (Gelci,  Chavz  and  Devillaz,  1963) 
and  the  Barnet  model  (Barnett,  1968)  are  also  being  tested 
and  evaluated  to  determine  which  model  would  provide  the  best 
results.  It  may  be  that  a  hybrid  model  containing  the  best 
features  of  each  model  may  evolve.  For  example,  the  Barnett 
model  has  a  wave-wave  interaction  term  which  may  improve  the 
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dissipation  of  wave  energy.  If  it  does,  such  a  term  will  be 
included  in  the  operational  model. 

The  initial  thrust  to  improve  the  present  operational 
model  will  be  slanted  towards  improving  the  wind  input.  Wind 
fields  will  be  derived  every  three  hours  rather  than  the  present 
six  hours  so  that  the  wind  input  can  coincide  with  the  basic 
time  step  of  the  wave  spectral  model.  Investigations  are 
being  made  to  determine  if  the  Cardone  Planetary  Boundary  Layer 
Model  modified  by  Dr.  Jack  Kaltala.?  FNWC  meteorologist-,  will 
produce  wind  fields  better  suited  for  use  as  input  to  the 
spectral  model  than  the  existing  FNWC  wind  analyses  and  prog¬ 
nosis.  It  is  believed  that  some  significant  improvement 
should  occur  since  the  wind  velocities  will  be  computed  for 
some  known  height  above  the  sea  surface  (preferable  19.5 
meters)  and  frictional  wind  velocities  (U*)  will  be  a  function 
of  atmospheric  stability. 

One  modification  that  would  probably  significantly  improve 
the  wave  spectral  model  is  to  increase  the  number  of  direction 
bands  used  to  define  the  directional  spectra  from  the  present 
twelve  bands  (30°  increments)  to  twenty-four  bands  (15°  in¬ 
crements)  .  This  would  better  define  the  travel  path  of  the 
wave  energy  over  long  distances.  Unfortunately,  doubling  the 
number  of  direction  bands  would  also  double  the  amount  of 
computer  storage  and  computation  time  needed  to  operate  the 
computer  model.  Thus,  for  operational  purposes,  this  modifi¬ 
cation  is  not  practical  on  the  present  FNWC  computers  and 
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must  be  held  in  abeyance  until  a  larger  computer  system  becomes 
available. 

A  similar  problem  exists  in  defining  the  frequency  spectra. 
It  has  been  suggested  that  for  certain  applications  of  the 
output  of  the  wave  spectral  model,  the  highest  frequency  should 
be  better  defined.  Essentially,  the  highest  frequency  band 
extends  from  0.164HZ  to  0.4HZ  which  could  be  subdivided  into 
two  or  three  smaller  bands;  however,  the  same  problem  as  was 
discussed  with  the  direction  bands  would  exist.  The  size  of 
the  two  dimensional  wave  energy  matrices  would  increase  and 
strain  the  present  FNWC  computer  resources.  One  short-term 
solution  would  be  to  delete  the  lowest  frequency  band 
(0.039HZ)  and  divide  the  highest  frequency  into  two  smaller 
bands.  This  would  not  increase  the  size  of  the  arrays  in  the 
computer  model.  Before  this  modification  can  be  made,  a 
determination  must  be  made  as  to  whether  the  deletion  of  the 
lowest  frequency  oand  will  significantly  alter  the  computer 
wave  spectra. 

The  NDBO  Bucy  wave  spectra  along  with  the  corresponding 
meteorological  data  will  be  closely  monitored  to  determine 
if  the  errors,  seen  as  described  in  Case  Studies  B  and  C  in  . 
the  computed  spectra  have  a  definite  trend  or  are  random. 
Measured  wave  spectra  from  other  sources  such  as  OSV  PAPA,  air¬ 
plane  lasers  and  satellites  will  also  be  monitored  for  verifi¬ 
cation  purposes. 
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In  addition  to  being  more  accurate  than  the  output  of  the 
singular  model ,  the  output  of  the  wave  spectral  model  has 
a  wider  range  of  applications  for  U.S.  Navy  use.  Investi¬ 
gations  are  now  underway  to  determine  the  feasibility  of 
applying  the  wave  spectra  to  such  projects  as  ship  motion 
studies,  acoustical  ambient  noise  problems,  and  shallow  water 
wave  refraction  studies. 

A  considerable  amount  of  effort  was  expended  by  a  number 
of  people-in  order  f or  the  SOWM  tro preach; the  status' If.  has 
today.  Consderably  more  effort  will  be  required  to  improve 
the  computer  model  and  to  make  full  use  of  its  potential. 
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APPENDIX  A 


EQUATIONS  USED  TO  GROW  AND  DISSIPATE  WAVE  ENERGY 
Comprehensive  derivations  of  most  of  the  following  equations 
can  be  found  in  Inoue  (1967)  and  Cardone  (1969). 

FNWC  analyzed  and  forecasted  surface  wind  speeds  are  used 
essentially  as  a  first  guess  in  the  following  equation  to  obtain 
frictional  wind  felocities  (U*)  and  surface  roughness  parameters 
(ZQ) .  Equation  (1)  is 


where  K  is  Von  Karman's  constant  and  has  a  value  of  0.4 
U  is  the  wind  velocity,  .  >  •  ■*-“* 

Z  is  the  elevation  at  which  the  wind  velocity  is 
measured  or,  in  this  case,  calculated, 

Zq  is  the  surface  roughness  parameter  and  in  this 
program  is  defined  as 

a/U*  +  +  T  and  a,  0  and  y  are  constants 

The  growth  equations  in  the  wave  spectral  model  are  designed 
to  use  6.1  meter  and  19.5  meter  wind  velocities.  If  the  winds 
input  is  based  on  some  other  height  than  those  above-mentioned, 
then  the  wind  speeds  should  be  recalculated  for  either  6.1  meters 
or  19.5  meters,  depending  on  which  equation  is  being  used.  The 
FNWC  wind  fields  are  only  considered  surface  winds  without  any 
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regard  to  height.  The  operational  wave  model  does  not  make 
any  attempt  to  correct  the  FNWC  winds  for  elevation. 

Since  there  are  two  unknowns  (U*,  ZQ)  in  equation  (1) ,  the 
equation  is  solved  by  reiterative  means.  Once  U*  and  ZQ  have 
been  calculated  equation  (1)  is  used  to  obtain  the  6.1  meter 
and  19.5  meter  wind  speeds. 

Neglecting  non-linear  effects/  the  spectral  components  can 
be  initially  expressed  as 


(2) 


where  A[f,U{t,x)  is  the  resonance  growth  mechanism, 

B(f,U(t,x)]  S{f,t,x)  is  the  Miles  instability  term, 

S (f , t, x)  is  the  spectral  density,  and 

f  is  the  wave  frequency,  t  is  time,  x  is  wave  direction 
Based  on  Cardor.e's  recent  work,  A(f,U)  can  be  defined  in  one 
of  two  ways  depending  on  the  magnitude  of  (27rf/U6,i).  If  the 
value  of  (2Trf/U6oi)  is  less  than  0.02,  then  A(F,U6#1)  is  defined 
as 


0.  675  *10~ 1 2 


fTl  2 


-IT  2 


0. 022  •  2  8  to4  U3 
_ §_*J- 

W  +  Y 


de  (3) 


where  W  -[0.270-:r2  +  {k  sin0}],  and. 

Y  =[V9  {0.02'  1*20}2,5+  {k  cos  (0-r) }  2  ] 
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Where  u>  ■  2irf, 


g  is  the  gravitational  acceleration, 
U  is  the  wind  speed  at  6.1  meters, 
k  =  u)2/g,  and 
r  =  w2/U6. i 


If  the  value  of  (2n/U6.i)  is  greater  than  0.02  then 
A(f,U)  is  defined  as 


'  TT2 


0.675.10 


1  2 


6.1 


l  ir2  0.2704  r2  +  IK  sin  6  ]2  [  V9  (-)  cos(6-r) }  2  ]  J 


d6 


f4) 


The  instability  growth  term  is  described  as  follows: 


B/F  = 


*r  {V.*k* 


p„  c^k  1  cos^e 

0) 


U" 

U’ 

I  [U  Cos  a-C]  £"kZdZ 

h  d 

Zm 

Z 

m 

(5) 


+  Ap  |ar2  (-U")  cos  6  ju  cos  0-C I  t“2kZdZ} 


Where  p  is  the  density  of  water, 

U) 

U"  is  the  mean  wind  profile  curvature, 

U'  is  the  mean  wind  profile  slope, 

8  is  the  directional  difference  between  the  wind  and  wave, 
A  is  a  constant  and  equal  to  tt, 
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A  is  a  constant  determined  by  experiment  and  ascer- 
P 

_  2 

tained  by  Phillips  to  be  1.6x10  with  an  uncertainty 
of  30%,  T2  is  +1  above  the  matched  layer  and  less 
than  +1  below  the  matched  layer. 

The  first,  term  on  the  right  hand  side  is  Miles'  solution 
and  the  second  term  is  a  later  contribution  by  Phillips.  B/f 
is  a  dimensionaless  quantity.  A  relationship  exists  between 
B/f  and  U*/C  (where  C  is  the  phase  velocity  of  the  wave  energy) 
and  a  number  of  investigators  have  examined  this  relationship 
on  a  theoretical  basis  or  in  field  experiments.  Inoue’s 
solution  which  is  used  in  the  wave  spectral  model  is 

B(f,U*)  =  {1.39  *  l(fV7“><>[U*/C  -  3.1  *  10~212  +  (6) 

+  0.725(U*/C) 2i"4  *  10"4 (C/U*) 2}f 

If  the  v;ave  growth  equation  was  used  in  the  computer  model 
without  any  constraints,  the  wave  energy  could  have  unlimited 
growth.  Since  this  does  not  occur  in  nature,  the  concept  of  a 
fully  developed  sea  was  introduced  in  order  to  approximate  the 
real  world.  The  essence  of  this  concept  is  that  if  a  wind  with 
the  same  magnitude  blows  in  the  same  direction  over  a  given 
fetch  for  enough  time,  the  wave  spectrum  will  become  fully 
developed  and  no  matter  how  much  longer  the  same  wind  velocity 
exists,  the  spectrun  will  no  longer  continue  to  grow.  The 
Pierson/Moskowitz  (P/M)  fully  developed  spectrun  is  used  to 
limit  wave  growth  and  has  the  form 
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smM 


(7) 


*  22I  e-e<“0/«)  “ 

w5 

where  a  =  8.1*10-3 
8  «  0.74,  and 
ioQ  =  g/U 

The  wave  spectrum  at  each  grid  point  is  compared  to  the  P/M 
spectrun  for  the  given  wind  speed.  The  P/M  spectrun  is  inserted 
in  equation  (2)  by 

[A  (1  •  f~  )2)1/2  +  BS]  [1  -  <|- )  2]  (8) 

*'ae  ^oo 

-T-  ■  -  ■  “  -  •  -  T  -  ' 

Ther  solution  for  zero  initial  conditions  is 
S (f , t)  =  A{&Bt-r,~  -1  >  [1  +  >2]“,/2  (9) 


Since  equations  (8)  and  (9)  only  provide  one  dimensional 
spectra,  an  equation  developed  by  the  Stereo-Wave  Observation 
Project  (SWOP)  is  used  to  obtain  the  directional  spectra.  This 
equation  is 


F  (w,  0  ,  U)  =  r(l*(0.5  =  0.82  C  '/z  (wU/g)  4  cos  2  6  (10) 


+  9. 32l~1/2  (wU/g)  4  cos  40j 
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for  -tt/2  <  6  <  tt/2  and  0  is  the  angle  between  the  wind  and 
wave  directions.  F  (to,  0  ,U)  is  equal  to  0  everywhere  else. 

The  directional  spectra  are  computed  for  30°  increments. 

As  an  example,  if  the  wind  direction  were  180°,  the  distribution 
of  wave  energy  would  be  as  follows:  37.5%  of  the  energy  would 
be  placed  in  the  180°  direction  band;  25%  would  be  placed  in 
the  150°  and  210°  bands;  and  6.25%  of  the  energy  would  be  placed 
in  the  120°  and  240°  direction  bands.  If  the  wind  direction 
were  190°  rather  than  180°,  then  the  energy  distribution  would 
be  more  biased  in  the  210°  and  240°  direction  bands,  than  in  the 
120°  and  150°  bands. 

Wave  energy  is  only  dissipated  when  the  waves  encounter 
land  or  interact  with  the  wind.  Inoue  developed  a  formula  to 
account  for  the  wave  dissipation  if  the  angle  between  the  wind 
and  wave  directions  is  greater  than  75°.  The  formula  is: 


where  is  spectral  component  after  dissipation, 

So(fi,M  is  spectral  component  before  dissipation, 
ff,©!  is  the  center  frequency  and  direction  of  the 
component, 

C  is  a  constant  690  for  ft2/sec  and  169.2  for  m2/sec) , 
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and  K < 0 1 )  is  0  for  6i  <  75° 

K (0 i )  is  1.5  for  75®  <  0  <  105° 

K (0 i )  is  3.0  for  105°  <  0j  <  135° 

K (0 i )  is  4.5  for  135°  <  0|  <  165°,  and 

K(0i)  is  6.0  for  165°  <  0j  <  180®. 

Significant  wave  height  (H^)  which  is  defined  as  the 
average  height  of  the  highest  one-third  waves  is  obtained  by 

5V,  =  4/E^  (12) 


where  E^ot  is  the  total  two-dimensional  wave  energy  at 
a  grid  point. 
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TABLE  1 


Central  Central  Frequency 

Frequency  (H2)  Period  (Seconds)  Bandwidth  (H  ) 


0.164 

6.1 

.164 

- 

to 

0.153 

6.5 

.142 

- 

.164 

0.133 

7.5 

.125 

~ 

.142 

0.117 

8.6 

.108 

-■ 

.125 

0.103 

9.7' 

.097 

- 

.108 

0.092 

10.9 

.086 

/ 

- 

.057 

li.  083 

12.0 

.080 

- 

.086 

0.078 

12.9 

.075 

.080 

0.072 

13.8 

.  .069 

- 

1.075 

0.067 

15.0 

.064 

- 

.069 

0.061 

16.4 

.058 

- 

.064 

0.056 

18.0 

.053 

- 

.058 

0.050 

20.0 

..047 

- 

.053 

0.044 

22.5 

.042 

- 

.047 

0.039 

25.7 

.036 

— 

.042 
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TA3LE  2 


PROPERTIES  OF  THE  GHOXONIC  PROJECTION  ON  A  FACE 
Area,  l/20th  of  the  Earth's  surface 
On  the  Earth: 

Length  of  side 

or 

Length  of  altitude 

or 

Vertex  angle  (spherical) 

On  the  Plane: 

Length  of  side 
Length  of  altitude 
Vertex  angle 


OF  AN  ICOSAHEDRON 
7.45  3C  10^  sq  n 

63°26.1' 

3806*  n  mi 

58°16.9' 

3497*  n  mi 
720 


4552.  n  mi 
3942.  n  mi 
60° 


Distortion  relative  to  1,00  at,  the  tangent  point: 

Radial  -  iwTrimm  (at  vertices)  1,58 

“  midpoint  of  sides  1,15 

Transverse  -  max?. run.  1.26 


-  midpoint  of  sides 
Areal  -  rax?  rr.;m. 

-  midpoint  of  sides 

Distortion  relative  tc  1.00  at  the  location  of 
mean  distortion: 

Radial  -  na.xl mum  (at  vertices) 

-  alnimna  (at  tangent  point) 
Transverse  —  maximum 

-  minima 
Areal  -  rax?,  mum 

-  minimum 


1.07 

1-99 

1.23 


1.40 

.88 

1.19 

.94 

1.65 

.82 


"Assumes  spherical  Earth  of  radius  3440.19  n  mi 


(Lockheed,  1966) 


mi 
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TABLE  3 


DTG 

YYMMDDHH 

LOCATION 

WIND  SPEED (K+S) 

WIND  DIRECTION 

SIGNIFICANT 
WAVE  HEIGHT (F: 

69112906 

48°N,165°E 

58 

320 

30 

69112912 

45°N, 170°E 

69-70 

260-270 

32-33 

69112918 

45°N, 170°E 

68 

290-300 

45 

6911300 

40°N,175°E 

72-73 

240 

50-52 

69113006 

45°N,178°E 

00 

vo 

1 

VO 

VO 

290-300 

t  ■ 

52 

69113012' 

:'40PN,176°W 

73-74 

250-260 

55-62 

69113018 

40~45°N,180°W 

87-91 

320-340 

74 

69120100 

40°N, 165°W 

78 

250-270 

c 

69-72 

69120106 

'42°H,170°W 

■  :  73  , ... 

330 

66-70 

„40PN,I65PW 

..  .  74- 

.  260 

69-75 
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TABLE  4 


WAVE  STATISTICS  AT  27°30'N,  157°45'WFOR  23Z,  1  DECEMBER  1969 


Largest  Wave  Height  in  One  Hour  Sample 
Significant  Wave  Height 
Average  Wave  Height 


24.4  meters 

14.4  meters 

10.5  meters 


80  feet 
49  feet 
34.5  feet 


Predominant  Wave  Spectral  Peak 


0.05-0. 062HZ  (20  to  16  sec) 


TABLE  5 


RATIO  OF  ACTUAL  WAVE  AMPLITUDE  TO  TUCKER  RECORDED  AMPLITUDE 


HERTZ  AMPLITUDE  RATIO 


.05 

.06 

.07 

.06 

.09 

.10 

.11 

.12 

.13 

.14 

.15 

.16 

.17 

.18 

.19 

.20 

.21 


1.06 

1.04 

1.04 

1.06 

1.09 

1.13 
1.18 
1.24 
1.32 
1.41 
1.52 
1.64 
1.78 
1.94 

2.14 
2.35 
2.62 


.22 

.23 

.24 

.25 

.26 

.27 

.28 

.29 

.30 

.31 

.32 

.33 

.34 

.35 

.36 

.37 

.38 

.39 

.40 

.41 

.42 

.43 

.44 

.45 

.46 

.47 

.48 

.49 

.50 


2.92 

3.28 

3.70 

4.20 

4.78 

5.48 

6.31 

7.31 
8.51 
9.96 

11.71 

13.85 

16.46 

19.66 

23.61 

26.49 

34.57 

42.16 

51.69 

63.69 
78.89 
98.23 

122.9 
154.7 

195.6 

243.6 

317.6 

407.9 
526.5 
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TABLE  6 


WIND  VELOCITY 


DTG 

GRD  PT  164 

GRD  PT  165 

OSV  PAPA 

50.9°N,145.6W 

48.5N.147.5W 

50°N,145°W 

KTS  DIR 

KTS  DIR 

73102512 

20.5,  227° 

22.5,  220 

73102518 

18.5,  216° 

19.5,  208 

73102600 

14.2,  207 

16.7,  206 

+3 

73102606 

15.5,  179 

20.1*  .200 

13,140115,1101 

73102612 

9.6,  142 

16.9,  213 

17,080117,090] 

73102618 

13.7,  201 

26.0,  219 

27,230[35,240] 

73102700 

25.0,.  228 

33.0*  229 

31 *210 [31 *200} 

73102706 

28.0,  209 

35.7/210 

31, 200 [38, 210] 

73102712 

25.9,  232 

30.6,  226 

35,220[32,200] 

73102718 

21.0,  233 

31.6,  224 

31,230[36,280] 

73102300 

25,0,  280 

30.8,  259 

36,310[30,300] 

73102806 

27.2  283 

32.6,  273 

27,280[I4,230] 

73102812 

21.8,  240 

26.2,  255 

24, 220[26,190] 

73102818 

34.3,  179 

34.7,  199 

37,180141,180] 

73102900 

35.6,  207 

38.9,  214 

34,200 
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TABLE  7 


COMPARISON  OF  SHIP  VISUAL  OBSERVATIONS  WITH  WAVE  SPECTRAL  MODEL 
SIGNIFICANT  HEIGHTS  AND  DIRECTIONS 


Ship 

DTG 

Model 

DTG 

LAT 

LONG 

Ship 

Speed 

Observed 

Height/ 

Direction 

SOWM 

Height/ 

Directior 

7502020UZ 

75020200Z 

37°  31 5  N 

124°38'W 

26 

13f t/270° 

15f t/305° 

75020216Z 

75020212Z 

37°15'N 

133°17'W 

26 

15f t/310° 

12f t/299° 

75020218Z 

75020212Z 

37° 12  *N 

134°30 '  W 

26 

17ft/340° 

16ft/305° 

75Q2022QZ 

75020300Z 

37° 12 ' N 

135°  30 ' W 

25 

20f t/350° 

12f t/305° 

75020223Z 

75020300Z 

37°12 'N 

137°36'N 

26 

13f t/320° 

12f t/321° 

7  5020302Z 

75020300Z 

37° 10  *N 

139°00*W 

26 

16f t/350° 

12f t/360° 

75020317  Z 

75020312Z 

36°  56 'N 

147°11'W 

27 

12f t/40° 
(200°) 

6ft/170° 

(15°) 

750203I9Z 

75Q20400Z 

36°  55 1 N 

148°20'W 

26 

lOf t/20° 
(190°) 

lOf t/172° 
(15°) 

75020321Z 

75020400 Z 

3  6  °  5  3  '  N 

149°20’W 

27 

8ft/20° 

(190°) 

lOf t/172° 
(15p ) 

75020323Z 

750204 00Z 

36°5,N 

150°24  ’W 

27 

8f t/220° 

lOf t/172° 
(15°) 

75020400Z 

7502040CZ 

36°  50  1 N 

150°56'W 

26 

8ft/210° 

lOf t/172° 

75020401Z 

75020400Z 

36°  49  1 N 

151°30'W 

26 

8ft/210° 

9f t/170° 

75020418Z 

7  502050CZ 

36°  34  '  N 

160°40'W 

21 

18f t/310° 

15f t/307  ° 

75020503Z 

750205G0Z 

36  °  22  1 N 

165°20 1 W 

26 

6ft/340° 

(290°) 

lOf t/300° 
(320°; 

75020519 

75020600Z 

36°05'N 

174°W 

26 

8f t/280® 

6ft/290° 

75020601Z 

75020600Z 

36°08  *N 

177°W 

26 

12f t/295° 

9ft/290° 

75020603Z 

75020600Z 

36°07'N 

178°09'W 

26 

12f t/290° 

9t t/290° 
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TABLE  7  -  CONTINUED 


Ship 

DTG 

Model 

DTG 

LAT 

LONG 

Ship 

Speed 

Observed 

Height/ 

Direction 

SOWM 

Height/ 

Direction 

75020622Z 

75020700Z 

35°  40 ' N 

172  °00 1 E 

26 

35ft/267° 

(270°) 

42f t/280° 

75020623Z 

75020700Z 

35°  40 ' N 

1 7  2  °  0  0  1 E 

9 

4 Oft/265 ° 
(27G° ) 

42f t/280° 

75020700Z 

75020700Z 

35°  36 ' N 

171°50'E 

3 

50f t/270° 
(280°) 

42f t/280  5 

75020702Z 

75020700Z 

35  0  36 ' N 

171°  35  1 E 

8 

30ft/275° 

(280°) 

40ft/280° 

75020704Z 

75020700Z 

35°  32 ' N 

171°13 *E 

8 

40ft/310° 

(300°')' 

40ft/?80° 
(270  6 ) 

7  5  0  2  0  7  0  5  Z 

75020700Z 

35 0  30  '  N 

171°03'E 

8 

40f t/290° 

40f t/280° 

75020721Z 

75020800Z 

3  4  °  4  0  '  N 

165°58 'E 

28 

25f t/350  0 
(310°)- 

2if t/330®  : 

7502072 3Z 

75Q20800Z 

34°30  '  N 

165°00 ' E 

31 

:  25f t/10°  : 

(190°) 

21ft/ 3*30"°"' 
(200°) 

75020801Z 

75C20800Z 

34 ° 30  ' N 

163°45’E 

31 

20ft/10° 

18f t/20  ° 

75020305 

75Q20300Z 

3  4  °  3  4  '  N 

161 °22 ' E 

30 

18ft/220° 

16ft/20° 

(200°) 

75020P06 

75020312Z 

34 7 3 4'  N 

160°44'E 

31 

25f t/200° 

18f t/10° 
(180°) 

75020822 

750209G0Z 

3 5°  32  1 N 

151°40 ' E 

30 

20f t/130° 

15ft/200° 

'  5  0  2  0  9  0  0  Z 

75020930Z 

3  5  0  2  5  ’  N 

150°20 ' E 

32 

14f t/180° 

15ft/200° 

7 5 020902 Z 

75020900Z 

3  5  5  2  5  1  N 

149° 10 ' E 

30 

lOf t/150° 
(180°) 

12f t/200° 

/5020907Z 

75020912Z 

35°10'N 

146°06 ' E 

30 

10ft/320o 

9f t/270  ° 
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ENERGY  GROWTH  (FIG1) 


ENERGY  GROWTH  (FIG  2) 


ENERGY  GROWTH  (FIG  3) 


d 
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PERIOD  (SEC)  IS. OO  9.00  8.4  3  5.00  4.09  3.45  3.00 

FREQ  (cyc/sec)  0.069  o.m  o.iss  0.200  0.244  0.299  0.33 


WAVE  ENERGY  GROWTH  (FIG  4) 
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WAVE  ENERGY  GROWTH  (Fig  7.) 
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USE  DIRECTLY. 


PRIMARY 


USE  LOCATIONS  AND  2  ALTERNATELY 
TO  GET  EQUIVALENT  OF  PATH  FROM 
LOCATION  3. 


SECONDARY 


FIG.  10 


Directions  to  Nearest  Neighbors  (Lockheed,  1966) 
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FIG.  II 

Bridging  Subprojections.  The  two  subprojections  are  shown  separated  along 
thoir  common  border.  The  incoming  energy  at  E  in  the  direction  shown  is 
propagation  from  the  interpolated  point  C  and  refracted  as  it  crosses  the 
border.  D  and  E  represent  the  same  points  but  have  their  energy  stored  in 
the  coordinates  of  subprojection  A-A  and  B-B  respectively.  (Lockheed,  1067) 
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WIND  SPEED  (m  /sec) 

FIG.  12 

Theoretical  wind  profiles  in  the  marine  surface  boundary  layer  for  a 

surface  stress  of  1  dyne/cm*  and  neutral  (N) ,  unstable  (u) ,  and  stable  (S) 
stratification.  (Cardone,  1969) 
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FIG.  29 

sowm  Two  Dimensional  Way#  Spectra 
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October  26  1973  J.8Z 


WAVE  ENERGY  GROWTH  (FIG  31) 


Ill  fmmr  in-. 


r  73  18  Z  WIND  FIELD 


FIGURE  34 

28  OCT  73  12  Z  WIND  FIELD 


WAVE  ENERGY  GROWTH  (FIG  37) 


EB03  12-24-74  OOOOZ 
56.  ON  I47.9W 
WIND*  282*  17.9  M/S 
H  1/3  »  78  M 
WAVE  PERIODS*  *  15  SEC 


FNWC  SPECTRAL  WAVE  MODEL 
12-24-74  OOOOZ 
56.24N  147. 54 W 
WIND- 263*  2QO  M/S 
H  I/3-S3M 

WAVE  PERI008<  15.0  SEC,  226* 


FIGURE  38 

(After  ODSI ,  1975) 


98 


EB03  12-26-74  OOOOZ 
56.  ON  147. 9 W 
WIND  331*  6.0  M/S 

HI/3-  3.1  M 
WAVE  PERIODS'  10.0  SEC 


FNWC  SPECTRAL  WAVE  MODEL 
12-26-74  OOOOZ 
56.24N  14754  W 
WIND  287*  9.6  M/S 
HI/3-  4.1  M 

WAVE  PER  100 S'  109  SEC,  258 


ES03  12-29-74  I200Z 
5a  ON  1479 W 
WIND  540*  96  M/S 
H  1/3  -  6.8  M 
WAVE  PERIODS 1  *  15  SEC 


PNWC  SPECTRAL  WAVE  MODEL 
12-29-74  1200 Z 

5624N  147  54 W 
WIND-  287*  IIJ  M/S 
H  1/3-  72  M 

WAVE  PERIODS'  13.0  SEC,  228* 


SHIP  ROUTE 

February  2  1975  00Z  -  February  9  1975  07Z 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWTlwi  Data  tntarad) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


recipient's  catalog  number 


mmx mwwA 


4.  TITLE  (and  Submit) 

An  Evaluation  o  :  a  Hemispheric  Opera¬ 
tional  Wave  Spectral  Model 

S.  TYRE  OF  REPORT  *  PERIOD  COVERED 

Final  Report 

S.  PERFORMING  org.  report  number 

7.  AuTHOAf*; 

Sheldon  M.  Lazanoff 

Norman  M.  Stevenson 

S.  CONTRACT  OR  GRANT  NUMBERfa) 

s.  performing  organization  name  and  address 

Fleet  Numerical  Weather  Central 

Monterey,  California  93940 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

’  1.  CONTROLLING  OFFtCC  NAME  AND  ADDRESS 

Fleet  Numerical  Weather  Central 

12.  REPORT  DATE 

July  1975 

Monterey,  California  93940 

12.  NUMBER  OF  PAGES 

103 

14  MONITORING  agency  N  AnE  »  ADDRESS^//  tlltltrtnl  from  Controlling  Oltlco) 

IS.  SECURITY  CLASS,  (ol  tbit  ttpon) 

Unclassified 

IS*.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 

'»  DISTRIBUTION  STATEMENT  tot  f hit  Rtport) 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  ■a:  rA»  tbttroel  tnttrt*  In  Blot*  30,  II  Blllaront  Itom  K opart) 

»•  supplementary  notes 

It.  KEY  *0*0$  fContlnuo  on  for—  ••  »W#  tt  n*c**«  *rv  • n4  fjmnttfy  by  block  mmbmr) 

spectral  ocean  wave  model  wave  growth 

wave  spectra  wave  propagation 

icosahedral  gnomonic  model  wave  dissipation 

wind  velocity  forecasting 

f  i c t i ona  1  wind,  v 


20.  ABSTRACT  (Contlnuo  on  rover* •  •  Ido  It  noeoooory  mt  Itontity  by  klbok  number) 

Through  the  cooperative  efforts  of  the  Naval  Oceanographic  Office 
and  Fleet  Numerical  Weather  Central  personnel,  the  open  ocean 
Icosahedral-Gnomonic  wave  spectral  computer  model  developed  by*  - 
Professes  Willard  Pierson  has  been  placed  in  real  time  operational 
use  at  FNWC  under  the  title  Spectral  Ocean  Wave  Model  (SOWM) .  The 
model  computes  semi-daily  on-time  analyses  and  forecasts  to  seventy 
two  hours.  A  gnomonic  projection  allows  great  circles  tc  be 


DD  1473  EDITION  OR  I  NOV  ••  It  OBSOLETE  UNCLASSIFIED  +&+A*- 

f/N  010 2* 014* (to I 


UNCLASSIFIED 


•BCIIBITV  Cl  ABBlPICATKM  A 


xrrnill] 


